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AFIT/GM/ENP/00M-06 

Abstract 

A wet-microburst event on 16 August 1994 at the Kennedy Space Center’s 
Shuttle Landing Facility alerted forecasters from the 45 Weather Squadron (45WS), the 
provider of weather support to the Kennedy Space Center (KSC) and Cape Canaveral Air 
Station (CCAS), to the challenges of wet-microburst prediction. Although there was no 
operational impact, this event caused the 45WS to revise their severe thunderstorm 
forecasting procedures to specifically address microbursts, resulting in the locally 
developed Microburst-Day Potential Index (MDPI). MDPI provides a several-hour 
outlook of microburst potential based on the results of the Microburst and Severe 
Thunderstorm (MIST) project. The 45WS also conducted a preliminary evaluation of the 
Wind INDEX (WINDEX) for the KSC/CCAS microburst forecast problem. WINDEX 
provides an estimate of the maximum observed gust speed that can be expected should a 
microburst occur. This thesis presents an evaluation of MDPI and WINDEX based on 
microbursts identified by Sanger (1999) in his KSC/CCAS microburst climatology. A 
new index for assessing microburst potential is also introduced, incorporating both the 
MDPI and WINDEX parameters. Overall neither the MDPI nor the WINDEX performed 
particularly well in this application. The MDPI showed very little improvement over 
random guessing, and the WINDEX showed very little correlation to observed maximum 
microburst gust speed. The new microburst potential index outperformed MDPI in 
almost all categories. Further refinement of the new index is needed to make it a more 
useful forecasting tool. 


IX 


AN EVALUATION OF MICROBURST PREDICTION INDICES FOR THE 


KENNEDY SPACE CENTER AND CAPE CANAVERAL AIR STATION 

(KSC/CCAS) 

1. Introduction 

1.1 Background 

High winds are a significant hazard to space-launch operations, including 
movement of the space vehicle to the launch pad, the launch. Space Shuttle landing, and 
takeoff and landing of the Space Shuttle ferry aircraft. An unforecasted high wind event 
during any of these phases can be catastrophic, causing numerous casualties and severe 
damage to the space vehicle and launch pad. Less visible than the launch are the many 
support operations that can be adversely impacted by unforecasted strong wind events. 
Examples of such operations include fueling/defueling of the space vehicle, launch tower 
repair, erecting vehicles and payloads, space vehicle transportation, and flying operations. 

One of the most difficult high wind events to predict is the downburst. A 
downburst is a strong thunderstorm downdraft that produces a “starburst” outflow of 
damaging wind at or near the surface. Downbursts have been sub-classified into two 
types: the macroburst and the microburst. A macroburst is a large downburst with 
damaging outburst winds extending horizontally more than 4 km (2.5 miles). Damaging 
winds from a macroburst can last from 5 to 30 minutes, with winds as high as 60 ms' 1 
(116 knots). Microbursts, on the other hand, are small downbursts with damaging 
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outburst winds not exceeding 4 km (2.5 miles) in horizontal extent. Microbursts 
generally last less than 10 minutes, with intense microbursts producing winds as high as 
75 ms' 1 (146 knots). Microbursts are further classified into “wet-microburst” and “dry- 
microburst” depending on the amount of precipitation associated with the downdraft 
(Fujita 1985). The small temporal and spatial scales of the microburst make them 
difficult to predict and pose a significant threat to space launch and support operations. 

On 16 August 1994, the Kennedy Space Center (KSC) and Cape Canaveral Air 
Station (CCAS) experienced several “downrush” wind events. The strongest wind gust, 
which was in excess of 33.5 ms' 1 (65 knots), was recorded at the Shuttle Landing Facility 
and was later attributed to a wet-microburst event (Wheeler and Roeder 1996). 

Forecasters assigned to the 45 th Weather Squadron (45 WS), the provider of weather 
support to Patrick Air Force Base (PAFB), the Kennedy Space Center, Cape Canaveral 
Air Station and the Eastern Range, forecasted airmass-type thunderstorm activity for that 
afternoon and expected high wind gusts but not near the magnitude that was experienced. 

Although there was no operational impact, this event sparked the 45WS to launch 
a full investigation into how such a severe wind event could elude forecasters. As a 
result of the investigation, the 45WS created a conceptual forecast model, the 
“Microburst Funnel,” to guide forecasters in microburst prediction (Figure 1). 

The Microburst Funnel begins with a review of thunderstorm and microburst 
climatology for the KSC and surrounding area. If forecasters predict thunderstorm 
activity in the next 6 to 10 hours microburst outlook techniques, such as the Microburst- 
Day Potential Index (MDPI), and the Wind Index (WINDEX), are used to determine if 
the environment is favorable for microburst development, and to determine the maximum 
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wind gust that can be expected should a microburst occur. Intermediate techniques, such 
as satellite imagery interpretation, are used to identify areas where thunderstorm 
development is imminent. Finally, nowcasting techniques, such as radar and visual 
identification are used just prior to issuance or non-issuance of a high wind warning. 



Figure 1 . Microburst Funnel (Adapted from Roeder 1999) 

This thesis focuses only on the two outlook techniques mentioned above, the 
MDPI and WINDEX indices, since most pre-launch operations, such as shuttle roll-out, 
can last up to eight hours. Early detection of a microburst producing environment in 
tandem with accurate microburst gust speed predictors can help mission planners 
postpone lengthy pre-launch operations when out-of-threshold winds are predicted, 
saving precious time and money. 
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Both MDPI and WINDEX are computed using upper-air sounding data. MDPI 
uses the difference between the minimum equivalent potential temperature (0 e ) found 
between 650 mb and 500 mb, and the maximum 0 e near the surface to assess the potential 
for microburst activity (Wheeler and Roeder 1996). WINDEX is an empirical index that 
gives an estimate of the maximum potential wind gust, at the surface, in knots (McCann 
1994). Both indices are useful since they give significant lead time (6-10 hours) and are 
easily computed using data that is readily available. 

Throughout this thesis, the terms “downburst” and “microburst” will be used 
interchangeably. This is common when the horizontal scale of the event is unknown 
(Caracena et al. 1990). All microbursts occurring at the KSC, CCAS, and Eastern Range 
will be assumed to be wet-microbursts since the thermodynamic requirements for dry- 
microbursts (an extensive, diy, sub-cloud layer) rarely exist on the Central Florida 
Atlantic Coast. 

1.2 Problem Statement 

Neither MDPI nor WINDEX have been evaluated by the 45WS for greater than a 
six month period. Additionally, neither index has been evaluated using the microburst 
cases identified by Sanger (1999) in his KSC/CCAS microburst climatology. Preliminary 
evaluation of MDPI indicated that A0 e values of greater than 30K suggest a high 
likelihood of a microburst, assuming thunderstorm or rainshower development. 

Evaluation of MDPI using a larger period of record is required. 

Preliminary evaluation of the WINDEX caused the 45WS to question its 
usefulness at its locale. Expansion of the period of record is required before WINDEX 
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can be ruled out as a gust speed predictor for microbursts along the Central Florida 
Atlantic Coast. The 45WS is also interested in synergy between these two indices. 

1.3 Objective 

This thesis seeks to evaluate the MDPI by extending the evaluated data set from 6 
months to 25 months (May-September 1994-1998). 

A Weibull distribution is fit to 282 microburst gust speeds identified by Sanger 
(1999) in his KSC microburst climatology (see Chapter 2, section 2.6). The Weibull 
distribution is integrated for each of the 45WS warning thresholds, providing a “first 
guess” probability of needing to issue a convective wind warning should a microburst 
occur. 

The WINDEX is also evaluated for its usefulness as a predictor of microburst gust 
speed for the Central Florida Atlantic Coast. Finally, a new index for assessing 
microburst potential is introduced. This new index incorporates both the MDPI and 
WINDEX variables. 

1.4 Importance of Research 

The 45WS provides weather support to launch and support activities conducted 
by the National Aeronautics and Space Administration (NASA), the Department of 
Defense (DOD), and other users of the CCAS, KSC, and Eastern Range launch facilities. 
Part of this support is issuance of convective wind warnings for the criteria listed in Table 
1. The desired lead-times listed in Table 1 are needed to ensure PAFB, KSC, CCAS, and 
Eastern Range personnel have adequate time to protect valuable resources from adverse 
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weather. Failure to meet these lead-times may result in property loss and possible loss of 
human life. Accurate wind gust prediction indices play a vital role in helping to meet 
these stringent requirements. 


Location 

Criterion 

Lead-time 


>35 knots 

30 minutes 

KSC (Sfc-300ft) 

>50 knots 

60 minutes 


> 60 knots 

60 minutes 

CCAS (Sfc-200ft) 

>35 knots 

30 minutes 


>50 knots 

60 minutes 

PAFB (Sfc) 

> 35 knots 

60 minutes 


>50 knots 

60 minutes 


Table 1 . Convective Wind Warning Criteria (Roeder 1999) 


1.5 Overall Approach 

Six steps were required to complete this research. First a Weibull distribution was 
fit to 282 microburst gust speeds identified by Sanger (1999) in his KSC microburst 
climatology. This distribution was integrated for each of the 45WS convective wind 
warning thresholds, providing a “first guess” of the probability of breaking a threshold 
should a microburst occur. Second was the identification of microburst-days and the 
identification of thunderstorm-days without microbursts. This was accomplished by 
identifying all thunderstorm-days for the KSC and CCAS between May - September 
1994-98 by analyzing surface observations from the KSC and CCAS. A day was 
considered a “thunderstorm-day” if the KSC or CCAS surface observation reported 
thunderstorms, rainshowers, cumulonimbus, or lightning in either the significant weather 
section or the remarks section. These days were cross-referenced with the KSC/CCAS 
microburst-days identified by Sanger (1999) in his KSC microburst climatology. Since 
Sanger did not identify microburst events for 1994 a computer program was written to 
identify microbursts occurring in May-September 1994. 
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Step three consisted of computing the MDPI and WINDEX parameters. CCAS 
sounding data furnished by the Air Force Combat Climatology Center was used to 
compute MDPI and WINDEX values. 

Step four consisted of a performance evaluation of both the MDPI and WINDEX 
equations. Statistical tools such as hit-rate, false-alarm rate, probability of detection, 
Heidke skill score, and the Pearson correlation coefficient (r) were used to measure the 
performance of each index. 

Step five was development of a new microburst potential index using the MDPI 
and WINDEX parameters by way of discriminant analysis. Lastly, a performance 
comparison between the original MDPI and the new microburst potential index was 
completed using the same performance measurements listed in step three. 

1.6 Organizational Overview 

Chapter 2 presents an overview of literature regarding convective downbursts, 
including a brief history of terminology and the development of the MDPI and WINDEX 
indices. Chapter 3 describes the data set used in this thesis and the computation of the 
MDPI and WINDEX parameters. Chapter 4 discusses the research methodology. The 
statistical techniques used for data analysis and development of a new microburst index 
are also covered. Chapter 5 summarizes the results and conclusions, and discusses 
opportunities for future research. 
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2. Literature Review 


2.1 A Brief History of Microburst Terminology 

How have meteorologists gone from downdraft to “microburst?” The best way to 
answer this question is to look at the events, or discoveries, that led to this term. The 
existence of downdrafts, associated with thunderstorms, has been known to 
meteorologists since the late 19 th century (Fujita and Caracena 1977). The most notable 
study of horizontal and vertical air currents in and around thunderstorms, the 
Thunderstorm Project, was conducted by Byers and Braham (1949) in the late 1940’s. 
Based on this project, Byers and Braham (1949) established the three stages of 
thunderstorm cells that are commonly taught at universities around the world: a) the 
cumulus stage characterized by updrafts throughout the cell, b) the mature stage with 
both updrafts and downdrafts, and c) the dissipating stage dominated by strong 
downdrafts. Downbursts are a special type of downdraft first categorized in the mid- 
1970’s. 

While investigating the Eastern Airlines Flight 66 crash that occurred at John F. 
Kennedy Airport, New York City on 24 June 1975, Dr. T. Theodore Fujita hypothesized 
that the accident was a result of a unique wind phenomenon that he called a “downbursf ’ 
(Fujita and Caracena 1977). Fujita defined a downburst as a downdraft with speeds 
“comparable to or greater than the approximate rate of descent or climb of a jet aircraft 
on the final approach or takeoff at 91 m (300 ft) above the surface.” He chose what he 
considered a conservative threshold value of 3.6 m s' 1 (12 ft s" 1 ) as the dividing line 
between an ordinary downdraft and a downburst. 
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In 1978, the first field research project on downbursts, the Northern Illinois 
Meteorological Research On Downburst (NIMROD), was conducted by the University of 
Chicago, and was headed by Fujita and Srivastava (Fujita 1985). Using three Doppler 
radars and 27 portable automated mesonet stations, the research team collected data on a 
large number of downbursts over a period of 42 days. Because of the overwhelming 
number of wind observations recorded during the 42 day period (approx. 1,632,960 
observations), Fujita (1985) developed a computer algorithm for identifying potential 
downbursts that occurred within the network (Table 2). Condition 1 sets a minimum gust 
speed of 10 ms' 1 (19.4 knots) to be considered a downburst. Conditions 2 and 3 state that 
the peak wind must be at least 5 ms' 1 (10 knots) faster than the five-minute mean wind 
speed prior to and after the peak wind. Conditions 4 and 5 ensure that the peak wind is 
not a gust “superimposed” on sustained high winds. Condition 6 excludes the gust front 
“which is characterized by an exponential decay of the gusty winds behind a front” 
(Fujita 1985). 


Condition 1 

Peak wind must be greater than 10 ms' 1 (19.4 knots) 

Condition 2 

Peak wind must be at least 5 ms' 1 (10 knots) faster than five-minute mean 
prior to peak 

Condition 3 

Peak wind must be at least 5 ms' 1 (10 knots) faster than five-minute mean 
after the peak 

Condition 4 

Peak wind must be 1.25 times the five-minute mean prior to the peak 

Condition 5 

Peak wind must be 1.25 times the five-minute mean after the peak 

Condition 6 

Five-minute mean prior to peak must be no more than 1.5 times the five- 
minute mean after the peak 


Table 2. Six Condition Algorithm for Identifying Potential Downbursts (Fujita 1985) 
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To obtain further evidence of a downburst, Fujita (1985) simultaneously plotted 
wind vectors for all towers at the time of the suspected downburst in order to analyze the 
field for a distinctive divergent flow pattern. A potential downburst was catagorized as a 
true downburst only if a divergent flow pattern was recognized. Fifty downbursts were 
identified during this project. 

Noticing significant differences in the horizontal scale of these downbursts, Fujita 
(1985) subclassified them into “macroburst” and “microburst.” A macroburst is a large 
downburst with damaging outburst winds extending horizontally more than 4 km (2.5 
miles). Damaging winds from a macroburst can last from 5 to 30 minutes, with winds as 
high as 60 ms" 1 (116 knots). Microbursts, on the other hand, are small downbursts with 
damaging outburst winds not exceeding 4 km (2.5 miles) in horizontal extent. 
Microbursts generally last less than 10 minutes, with intense microbursts producing 
winds as high as 75 ms" 1 (146 knots) (Fujita 1985). 

Of the 50 microbursts detected during the NIMROD field project, 64% were 
accompanied by 0.25 mm (0.01 inches) of rain or more. The remaining 36% had no 
measurable rain. Because of this, microbursts were further classified into “dry 
microbursf ’ and “wet-microbursf ’ depending on the accompaniment of measurable rain 
(Fujita 1985). Since then, research by Brown et al. (1982), Caracena et al. (1983), 
Wakimoto (1985), and many others, has helped to increase understanding of the dry 
microburst environment. Wet-microbursts have not been studied nearly as much; 
however, it is clear that the wet-microburst environment is quite different from the dry- 
microburst environment (Caracena et al. 1990). 
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2.2 Wet-Microburst Storm Structure 


In 1975, a meteorological mesonetwork was established in south central Florida 
for the Florida Area Cumulus Experiment (FACE). The FACE project was not a 
microburst study; however, on 1 July 1975, a wet-microburst made nearly a direct hit on 
the experiment’s Field Observing Site (FOS). Weather observers at the FOS originally 
believed that a tornado had touched down nearby. An investigation shortly after the 
storm determined that the damage pattern was indicative of outflow at the base of a 
downdraft versus flow into a tornado vortex (Caracena and Maier 1987). Because of the 
fine data resolution of the collected data, and the central location of the microburst within 
the mesonetwork, this microburst, known as the FACE microburst, has significantly 
enhanced current understanding of wet-microburst structure (Caracena and Maier 1987). 

Generally, thunderstorms producing wet-microbursts are characterized by low 
cloud bases, storm tops reaching as high as 15 km (49,000 ft) and strong precipitation 
cores that are composed almost entirely of ice in the upper-levels (Atkins and Wakimoto 
1991). They normally occur in environments that are nearly moist adiabatic and 
statically unstable from the surface to about 500 mb, with a deep layer of high relative 
humidity near the surface (Doswell 1994). This rather moist layer is capped by a cool, 
dry layer (low equivalent potential temperature) which is the most probable source of 
microburst energy once convection and precipitation begin to occur (Caracena and Maier 
1987). As downdrafts from precipitation loading begin, the cool, dry layer ejects pockets 
of dry air into the saturated downdraft region (Caracena et al. 1990). Latent heat of 
evaporation into the dry air contributes significantly to the negative buoyancy of the 
parcel (Caracena and Maier 1987). The key to developing a wet-microburst is this 
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evaporation of condensed water during descent, thus maintaining saturation. Without this 
evaporation, the parcel would begin to warm at the dry adiabatic lapse rate, become 
warmer than its surrounding environment, and quickly lose its negative buoyancy 
(Doswell 1994). Figure 1 shows a schematic of a thunderstorm producing a wet- 
microburst. In addition to the evaporative effects, the presence of high relative humidity 
in the low levels increases the virtual temperature difference between the downdraft and 
the environment, also contributing to an increase in negative buoyancy (Atkins and 
Wakimoto 1991). 


Wet-Microburst 



Figure 2. Conceptual model of a wet-microburst characterized by a dry source layer that 
ejects dry pockets of air into an underlying saturated layer, producing the evaporation that 
can cause a microburst (After Caracena et al. 1990). 
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2.3 The Wet-Microburst Environment 


Using data collected during the Microburst and Severe Thunderstorm (MIST) 
project that operated in northern Alabama in 1986, Atkins and Wakimoto (1991) 
documented the general environmental conditions (thermodynamic characteristics) that 
favor wet-microburst development. The MIST project used a network of 41 Portable 
Automated Mesonet (PAM-II) stations and 30 FAA-Lincoln Laboratory Operational 
Weather Studies (FLOWS) stations, spaced approximately 2 km (1.25 miles) apart, to 
collect meteorological data every minute. Atmospheric soundings were also taken twice 
daily, one at 0700 CDT, and one at 1300 CDT. From this set of data, 33 microbursts 
were identified using a computer algorithm first suggested by Fujita (1985). Analyzing 
the synoptic charts, sounding data, and mesonetwork data from the MIST network area, 
Atkins and Wakimoto (1991) determined the thermodynamic characteristics for 
microburst-producing days and compared them to the characteristics for the days with 
thunderstorms that did not produce microbursts. 

Atkins and Wakimoto described a wet-microburst environment as one that begins 
with a shallow morning radiation inversion that inhibits convection in the lowest 75 mb 
of the atmosphere. As the surface begins to heat, the inversion is replaced by a dry- 
adiabatic sub-cloud layer that extends from the surface to about 850 mb. The layer 
between 850 mb and 500 mb is relatively moist. This is capped by a cool, dry layer 
above 500 mb (Atkins and Wakimoto 1991). Figure 2 shows morning and afternoon 
sounding models for a wet-microburst-producing day. 
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Figure 3. Model Thermodynamic Charts Conducive for Wet-Microburst Occurrence. 
Dashed line is the dew point curve; solid line is the temperature curve (After Atkins and 
Wakimoto, 1991). 

By plotting the equivalent potential temperature (0 e ) profiles for the microburst 
producing days, Atkins and Wakimoto found a minimum of 0 e was typically found 
between 650 mb and 500 mb, or just below the capping cool, dry layer. The afternoon 0 e 
profile showed a A0 e , from the surface to the mid-level minimum, greater than or equal to 
20 K. The non-microburst producing storms had a A0 e less than or equal to 13 K. 
Equivalent potential temperature plots for other well-documented wet-microburst cases 
confirmed these results. Atkins and Wakimoto then suggested that by forecasting an 
afternoon sounding based on the morning sounding, the A0 e values of 20 K and 13 K can 
be used as “threshold” values for the development of wet-microbursts. They envisaged 
these thresholds being used to issue a general “area-wide alert” for potential wet- 
microburst activity. Figure 4 shows model morning and afternoon 0 e profiles for a wet- 
microburst environment. 


14 
















i i i n i 

| | | y j Morning 


i i i i / i 

||!l! Afternoon 

400mb J _L 1 / J _1 -L ]_ 

1 1 1 / 1 1 1 | ! 


_ 400mb j _[ j/L J _i x i_ 

1 1 1 / 1 1 III 

| | l/l 1 III 

i i y \ i i l f 

i i/i i ii i i 

500mb 1 l/l 1 1 1 1 1 


1 1 1/ 1 1 1 | 1 

1 1 / 1 1 1 1 1 

500mb 1 1 /I 1 1 III 

- t —- 1 -1-- f. — 1 - 1 - 

1 / 1 1 1 III 


- f - 1 —L i - 1 - y -*t--|- 1 - 

1 1 / 1 1 1 III 

i /I i i i iii 

_600mb_J_/|__[__|_l_L_-I_1_ 


1 l/l 1 1 1 | 1 

600mb 1 W 1 1 1 * 1 1 

IV i i i i ill 

__ i V i i i i i i i 

700mbj L L J JL _L L 

1 ixj 1 1 i 1 i 


i /T iii iii 

700mb j Vi l L J 1 L L 

1 Tv. Ill III 

_M0mb_|_[__l_i_1_ 

1 1 1 1 — . 1 1 1 1 


QAn . 1 1 1 1 1 1 1 
_800mb j_ i __l__|__J_ 

i i i i i i i 

HftA - U •III ^fx | ) 

_ ZMUUJLL—1_1_J_1_l-A 1-J-1- 


nnn«,u • 

^ U11,U _i- _L -L L _L xi. 4- L 

_!Qp0mbl _ 1 _| j _ | _ j_ 


lOOOrnb 1 1 1 1 1 0_L 1 



-j- -r r r -I -f r r 


Figure 4. Model 0 e profile (solid curve) for an environment conducive to wet-microburst 
development. Vertical dashed lines are 0 e values increasing from left to right at a 5 K 
interval (After Atkins and Wakimoto 1991). 

2.4 The Microburst-Dav Potential Index fMDPI) 

A wet-microburst event that occurred at the Kennedy Space Center’s (KSC) 
Shuttle Landing Facility on 16 August 1994 led the 45WS to revise their severe 
thunderstorm forecasting procedures to specifically address microbursts. This resulted in 
the development of the Microburst-Day Potential Index (MDPI). MDPI provides a 
several-hour outlook of microburst potential based on 0 e profiles first introduced by 
Atkins and Wakimoto (1991). The MDPI equation is: 

, _ T Maximum 6e - Minimum 6 e 

MDPI =- (lj 

CT 

where 

Maximum 0 e \s the maximum 0 e in the lowest 150 mb of the sounding (Kelvin) 

Minimum 6 e is the minimum 0 e in the layer between 650 and 500 mb (Kelvin), and 
CT is the critical threshold (or normalization factor) currently defined as 30 K. 
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Given thunderstorm development, MDPI values of 1 or greater indicate a strong 
probability of a wet-microburst, and MDPI values of less than 1 indicate a low 
probability of a wet-microburst. Preliminary analysis of MDPI indicated that it shows 
good skill in alerting forecasters to the potential for wet-microbursts (POD = 96.4%) 
without an unreasonable false alarm rate (FAR = 32.5%) (Wheeler and Roeder 1996). 

2.5 The Wind INDEX (WINDEX) 

McCann (1994) introduced an index for identifying airmasses favorable for the 
development of microbursts. The “Wind INDEX,” or WINDEX, is based on studies of 
observed and modeled microbursts. It is a highly empirical index based on the vertical 
momentum and continuity equations developed by Wolfson (1990). Wolfson showed 
that microburst vertical velocity is proportional to some “forcing” times the depth of the 
downdraft (Az). “This forcing is proportional to the square of the environmental lapse 
rate, T, through which a downdraft descends (McCann 1994).” Since a major source of 
microburst negative buoyancy is from the absorption of latent heat due to melting and 
evaporation, McCann considers the downdraft depth to be the height of the melting level 
above the surface. The WINDEX equation is: 

WI = 5[H m Rq(- r* - 30 + Q l - 2 Qm)] 0 5 (2) 

where 

Hm is the height of the melting level above the ground (km) 

Wis the lapse rate from the surface to the melting level (°C km' 1 ) 

Ql is the mixing ratio of the lowest 1 km above the surface (g kg' 1 ) 

Qm is the mixing ratio at the melting level (g kg' 1 ) 
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N is the moisture adjustment 12 g kg' 1 , and 
Rq = QiJN but not exceeding 1. 

Before using WINDEX, it is important to consider the following: 1) When the 
lapse rate is low, the radicand may become negative. When this happens, WI is zero and 
the probability of a microburst is nil. 2) The term Rq attempts to account for an 
overestimation of WINDEX values in a dry environment. This term implies that when 
the low-level mixing ratio is less than 12 g kg' 1 , the atmosphere is too dry to produce a 
high precipitation storm. 3) WINDEX is very sensitive to environmental lapse rate error. 
4) Multiplication of the radicand by five, on the right side of the WINDEX equation, 
gives an estimate of the maximum potential surface gust in knots (McCann 1994). 

WINDEX is computed using observed sounding data or data from numerical 
prediction models. The actual WINDEX value is the maximum potential surface gust, in 
knots, that can be expected from a microburst (McCann 1994). Based on McCann’s 
verification data, WINDEX correlates well with the observed gust speeds of known 
microburst events. Table 3 shows WINDEX values and observed gust speeds for well 
documented microbursts. 


LOCATION 

DATE 

WINDEX 

OBSERVED 

SOURCE 

Southern FL 

1 Jul 75 

53 

60 

Caracena and Maier (1987) 

Northern AL 

20 Jul 86 

59 

60 

Wakimoto and Bringi (1988) 

Northern AL 

13 Jul 86 

56 

56 

Wakimoto and Bringi (1988) 

San Antonio, TX 

1133191 

63 

59 

Ladd (1989) 

San Antonio, TX 


68 

53 

Ladd (1989) 

Norman, OK 


54 

55 

Stewart (1991) 

Tampa, FL 

19 Jul 88 

56 

52 

Stewart (1991) 

Vero Beach, FL 

4 Jun 90 

61 

50 

Stewart (1991) 

Amarillo, TX 

21 May 86 

44 

47 

Sohl (1987) 

near Denver, CO 

19 May 82 

47 

48 

Wakimoto (1985) 


Table 3. WINDEX Calculations for Known Microburst Events (McCann 1994). 
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Preliminary analysis of WINDEX by the 45WS showed that it tends to 
underestimate gust speeds for the Central Florida Atlantic Coast (Roeder 1999). Further 
analysis is required to determine the usefulness of WINDEX at the KSC and CCAS. 

2.5.1 Forcasting Dry Microbursts Using WINDEX 

Murdoch (1997) evaluated WINDEX as a gust predictor for dry-microbursts in 
the Southwest United States. Dry-microbursts were identified using Storm Data (1983- 
1993) and other dry-microburst studies. Nineteen upper-air soundings taken at 12 UTC 
between the months of March - October, from Midland, Amarillo, and El Paso, Texas and 
Albuquerque, New Mexico were loaded into the Skew-T/Hodograph Analysis and 
Research Program (SHARP) workstation for analysis. The 19 soundings were provided 
by the National Climatic Data Center (NCDC) and from the National Weather Service 
Office (NWSO) Midland’s upper-air data archives. Once loaded into SHARP, the 
soundings were modified to reflect the conditions at the time of the microburst. These 
“modified” soundings were used to calculate WINDEX. 

Murdoch (1997) concluded that WINDEX “agreed well with the observed and 
estimated wind gusts, yeilding a mean error of -5.8 [knots].” Murdoch (1997) offered 
two explanations for this underestimation in cases where the cloud bases were high: 1) 
Loss of momentum as the absorption of latent heat ends and adiabatic heating continues, 
and 2) Low average 1 km mixing ratios compared to a typical dry microburst sounding. 
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2.5.2 WINDEX Computation from Weather Satellite Derived Soundings 


The National Environmental Satellite Data and Information Service (NESDIS) 
Forecast Products Development Team (FPDT) evaluated two experimental microburst 
prediction products derived from Geostationary Operational Environmental Satellite 
(GOES) sounder retrievals (vertical temperature and moisture profiles) (Ellrod and 
Nelson 1998). One of the experimental products plots WINDEX as a color-coded 
graphic superimposed on a visible, infrared, or water vapor image. This product was 
evaluated for a two-month period. WINDEX values were compared to reports of wind 
damage or wind gusts greater than 50 knots reported in the National Weather Service 
(NWS) Storm Prediction Center’s (SPC) preliminary storm data. More than 300 wind 
damage reports were compared to adjacent WINDEX values derived from sounder 
retrievals. In 92 percent of the cases WINDEX values were greater than 40 knots, while 
60 percent exceeded 50 knots. Overall, Ellrod and Nelson (1998) found that the absolute 
value of the mean GOES WINDEX differed from the mean observed wind gusts by less 
than 3 knots. They found that the GOES WINDEX product was “generally most reliable 
over the eastern and central United States during daylight hours” (Ellrod and Nelson 
1998). 

2,6 Microburst Identification at the Kennedy Space Center (KSC) 

A four-year, summertime microburst climatology was produced for the KSC, 
CCAS, and Eastern Range by Sanger (1999). This is the first microburst climatology 
ever produced in the United States. The period of record was May through September 
1995-98. A total of 282 microbursts were identified over 114 microburst days during this 
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four-year period. Sanger showed that the most prominent months with microbursts, 
within the KSC WINDS network, are June, July, and August with the highest frequency 
in July. “The most favorable time for microbursts is between 1600 UTC (12 P.M. EDT) 
and 2200 UTC (6 P.M. EDT) with the peak occurring between 2000 UTC (4 P.M. EDT) 
and 2200 UTC (6 P.M. EDT)” (Sanger 1999). Appendix A lists the microburst events 
identified by Sanger. 

2.6.1 The KSC/CCAS Weather Information Network Display System (WINDS) 

The KSC operates one of the densest mesonetworks in continuous operation in the 
United States. This network collects, disseminates, and archives wind direction and 
speed, as well as other meteorological parameters, over an area of approximately 1200 
km which encompasses the KSC, CCAS, Eastern Range, and surrounding remote 
locations. On average there is one tower every 27 km 2 . The most densely instrumented 
region is along the East Coast within the KSC, CCAS, and Eastern Range complexes 
(Figure 5). The instrumented towers that make-up this network take meteorological 
measurements every minute, but archive data at five-minute intervals. Therefore, the 
data used for Sanger’s (1999) research consisted of the peak and mean wind speeds every 
five-minutes verses the one-minute interval used in the NIMROD and MIST studies 
(Sanger 1999). 

Sanger (1999) mentions several limitations of the WINDS network that make 
microburst identification difficult. First is the overall placement of the towers. Sanger 
mentions that many of the towers have “poor meteorological exposure” resulting in 
erroneous data. Second is the overall tower spacing. Towers to the south and southwest 
of CCAS are sparsely populated potentially inhibiting detection of a microburst between 
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towers. Third is the varying number of sensors on each tower. Some towers may only 
have sensors at 6 and 30 feet (Tower 512) while others may have as many as 8 sensors at 
6, 12, 54, 162, 204, 295, 394, and 492 feet (Tower 3131). Lastly, Sanger believes that the 
data interval of five-minutes may be too long for identifying short-lived microburst 
events (Sanger 1999). 

2.6.2 Microburst Identification Technique 

Sanger (1999) used a computer program, written in the Interactive Data Language 
(IDL), to analyze WINDS tower data from May through September 1995-98 to identify 
possible microbursts. Six conditions (Table 4) needed to be met simultaneously before a 
wind reading could be considered a possible microburst event. These conditions are a 
modified version of Fujita’s (1985) six condition algorithm for identifying potential 
downbursts. Since the WINDS archive data interval is 5 minutes, the pre-peak mean 
wind speed and post-peak mean wind speed used by Fujita were replaced by the peak 
wind 5 minutes before and 5 minutes after the potential microburst. 


Condition 1 

Potential microburst must be greater than 10 ms' 1 (19.4 knots) 

Condition 2 

Potential microburst must be at least 5 ms' (10 knots) faster than the five- 
minute peak prior to the potential microburst 

Condition 3 

Potential microburst must be at least 5 ms' 1 (10 knots) faster than the five- 
minute peak after the potential microburst 

Condition 4 

Potential microburst must be 1.25 times the five-minute peak prior to the 
potential microburst 

Condition 5 

Potential microburst must be 1.25 times the five-minute peak after the 
potential microburst 

Condition 6 

Five-minute peak prior to potential microburst must be no more than 1.5 
times the five-minute peak mean after the potential microburst 


Table 4. Sanger’s Six Condition Algorithm for Identifying Potential Microburst (Sanger 
1999) 
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Figure 5. Map of KSC WINDS Tower Network. Triangles indicate tower location. 
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Like Fujita, Sanger simultaneously plotted wind vectors for all towers at the time 
of the suspected microburst in order to analyze the field for a distinctive divergent flow 
pattern. Sanger, however, did not throw-out a suspected microburst if a divergent flow 
pattern could not be identified, noting that the wide spacing of towers could miss the 
divergent flow pattern. This is consistent with the technique used in the 1986 MIST 
project (Atkins and Wakimoto 1991). 
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3. Data and Computation of MDPI and WINDEX Parameters 


3.1 Introduction 

Three data-sets were used to accomplish this research: a) Weather Information 
Network Display System (WINDS) tower data from May - September 1994, b) Upper- 
air sounding data from May - September 1994-98, and c) KSC and CCAS surface 
observations from May - September 1994-98. The following gives a brief description of 
these data-sets and an overview of how they were used in the context of this research. 

3.2 Weather Information Network Display System (WINDS) Tower Data 

The Weather Information Network Display System (WINDS), operated by the 
KSC, collects, disseminates, and archives data from 201 meteorological sensors on 44 
instrumented towers at the CCAS and KSC, and at remote surrounding sites. This 
network is one of the most dense mesonetworks in continuous operation in the United 
States, covering approximately 1200 km for an average of one tower every 27 km . 
These towers are organized into three groups based on their primary application 
(Raytheon 1996): 

a. Launch Critical Towers - These towers are at the launch complexes, the 
Shuttle Landing Facility (SLF) and in some remote locations. They use the most 
sensitive sensors in the network and are used to ensure that launch constraints are 
satisfied during countdown and major pad operations. These sensors are also used to 
assess possible blast damage effects due to detonation of the solid rocket motors. 
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b. Safety Critical Towers - These towers are located near areas where propellants 
and other toxic chemicals are stored or handled. They are used by Range Safety models 
to predict the diffusion of potential airborne contaminants released from chemical spills. 

c. Forecast Critical Towers - These towers generally surround the KSC and 
CCAS, and are used primarily for weather forecasting and early detection of hazardous 
weather conditions. 

It is important to mention that, although each tower is assigned to one of the 
above groups, the data from these towers is displayed and archived as one integrated 
network, allowing any tower to contribute to any application. The distinction between 
tower groups is, therefore, insignificant operationally (Raytheon 1996). 

Of primary concern for this research are the wind speeds and gust speeds recorded 
by this network. WINDS data for May - September 1994 was provided by the 45WS in 
the format shown in Appendix B. The following parameters were provided: Julian Day, 
time (UTC), tower number, sensor height (ft), mean wind direction (degrees), 5-minute 
average wind speed (knots), instantaneous peak speed recorded in a 5-minute interval 
(knots), standard-deviation of wind direction (degrees), temperature (°C), dew-point 
temperature (°C), and relative humidity (%). Not all sensors report all of the above 
parameters. This data was used to identify potential microbursts that occurred in 1994, 
using Sanger’s six condition algorithm for identifying potential microbursts (Sanger 
1999). 

3.3 KSC Upper-Air Sounding Data 

Upper-Air Sounding Data from the KSC (station 74794) for the period May - 
September 1994-98 was provided by the Air Force Combat Climatology Center. 
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Soundings are taken routinely during the summer months at the KSC, at 10 UTC, 15 
UTC, and 22 UTC. Additional soundings are taken on launch days, and oh days with 
highly weather sensitive operations. The MDPI and WINDEX parameters were computed 
using the 15 UTC soundings since most microbursts, at KSC, occur between 16 UTC and 
22 UTC (Sanger 1999). The 10 UTC sounding is too early and would have little 
predictive ability for microbursts that occur after 16 UTC. Use of the 22 UTC sounding 
would have virtually no predictive ability since most microbursts occur before this time. 
The sounding data provided includes the following: date, time, geopotential height, 
pressure level, temperature, dew-point temperature, wind direction, and wind speed. 
Missing data is indicated by a “999” entry. Appendix C is a sample upper-air sounding in 
decoded format. Column 1 is the time (UTC); column 2 is the day; column 3 is the 
month; column 4 is the year; column 5 is the pressure (mb); column 6 is the height 
(meters); column 7 is the temperature (°C); column 8 is the dew-point temperature (°C); 
column 9 is the wind direction (degrees); column 10 is the wind speed (knots). 

3.3.1 Handling Missing Data Within a Sounding 

Many of the WINDEX and MDPI parameters require temperature and/or dew¬ 
point temperatures for levels that were reported as missing (reported as “999”). A 
computer program, written in IDL, was used to interpolate this data (Appendix D). 
Interpolation was done by assuming a constant lapse rate from the last known 
temperature (or dew-point) to the next known temperature (or dew-point) in the 
sounding. Figure 6 illustrates this interpolation process. 

Step 1 takes the difference between the last known and next known temperatures 
in the sounding. Step 2 computes the difference between the height of the last known 
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temperature value and the next known temperature value in the sounding. Step 3 builds a 
weighting ratio by taking the ratio of the difference between the height of the missing 
temperature value and the value computed in Step 2. This weighting ratio ensures that 
the temperature used in place of the missing value is proportional to the change in height 
from the last known temperature. Step 4 multiplies the weighting ratio by the value 
computed in Step 1 and subtracts that value from the last known temperature, yielding the 
interpolated temperature value. This same process is used to interpolate missing dew 
point temperatures as well. 




Figure 6. Interpolation of Missing Temperature Data. A weighting ratio was computed 
and then multiplied by the change in temperature between the two known levels. This 
value was then subtracted from the last known temperature reading. This ensured that the 
temperature was reduced by an amount that was proportional to the change in height from 
the last known temperature to the height of the missing value. 
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3.3.2 Computation of the MDPI parameters 

The Microburst-Day Potential Index uses the difference in equivalent potential 
temperature between the surface (lowest 150 mb) and mid-levels (between 650 mb and 
500 mb) to determine the potential for microburst activity within the next 6-10 hours. 
Equivalent potential temperature (0 e ) is the temperature a parcel of air would be if it was 
expanded pseudo-adiabatically until all the water vapor has condensed, released its latent 
heat, and fallen out, and was then compressed dry-adiabatically to 1000 mb (Wallace and 
Hobbs 1977). The following equations, from A WS/TR-83/001 Equations and Algorithms 
for Meteorological Applications in Air Weather Service (Duffield and Nastrom 1983), 
were used to compute equivalent potential temperature: 


= 71 


( u \*(l-0.28w)' 


J 


exp 


Tlcl 




-2.54 


M^l + O.Slw) 


J 


(3) 


where 


Tlcl = Td-(B + 0.00157 \Td - 0.000436^X71. - Td) + C 


(4) 


and 

d e is the equivalent potential temperature (Kelvin) 

T is the temperature at pressure level P (Kelvin) 

P 0 is the reference pressure level 1000 mb 
P is pressure (mb) 

A is the empirically derived value 3376 K 

Tlcl is the temperature at the lifted condensation level (LCL) (Kelvin) 
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A: is 0.2854 


w is the mixing ratio (kg kg' 1 ) 

B is the empirically derived value 0.212 °C 
C is the value 273.16 K which converts Celsius to Kelvin 
T c [ is the surface dew point temperature (°C), and 
T 0 is the surface temperature (°C) 

A computer program written in the Interactive Data Language (IDL) was used to 
compute 9 e at each level of the atmospheric sounding (Appendix D). The following 
equations were used to compute the mixing ratio (w) needed in equation (3): 

w= 0.622-—^— (5) 

p-e 


where 


e(Td ) = D • exp 


17.67-7* 
Td + E 


( 6 ) 


and 

w is the mixing ratio (kg kg' 1 ) 

p is the pressure at the level where T is measured (mb) 
e is the vapor pressure at the level where 7* is measured (mb) 

D is the empirically derived value 6.112 mb 
E is the empirically derived value 243.5 °C 

Td is the dew-point temperature at the level where T is measured (°C) 
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Equation (6) approximates the vapor pressure, and is valid within 0.1% over the 
temperature range -30°C <Td< 35°C (Rogers and Yau 1996). 

3.3.3 Computation of the WINDEX Parameters 

WINDEX uses five parameters that are computed from an atmospheric sounding, 
to estimate the maximum gust speed at the surface that can be expected from a 
microburst. These parameters are the height of the melting level (H m ) (km), the mean 
environmental lapse rate from the surface to the height of the melting level (/") (°C km' 

! ), the mixing ratio in the lowest 1 km of the atmosphere ( Qi ) (g kg' 1 ), the mixing ratio at 
the melting level (Q m ) (g kg' 1 ), and a moisture adjustment factor (R q ) which is equal to 
Q[/12 but not greater than 1. A computer program written in IDL was used to compute 
WINDEX (Appendix D). The following is a brief overview of how each parameter was 
computed: 

a. Height of the Melting Level (H m ) - The height where the temperature sounding 
first becomes less than zero (McCann 1999). This height was converted to kilometers by 
dividing by 1000. 

b. Mean Environmental Lapse Rate (/") from the Surface to H m - This was 
computed by dividing the surface temperature by the height of the melting level ( H m ) 
(McCann 1999). 

c. Mixing Ratio in the Lowest 1 km (Qi) - Computed by using equations (5) and 
(6), above, to compute the mixing ratio at the surface and at 1 km. The arithmetic mean 
was then computed and used as Qi. 

d. Mixing Ratio at the Melting Level ( Q m ) - Computed by using equations (5) and 
(6), above, at the level where the temperature sounding first becomes less than zero. 
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e. Moisture Adjustment Factor ( R q ) - This is equal to QJ 12 but not greater than 1. 
3.4 KSC and CCAS Surface Observations 

Surface observations from the KSC (KCOF) and the CCAS (KTTS), for the 
period May - September 1994-98, were used to identify thunderstorm days. A day was 
considered a “thunderstorm-day” if the KSC or CCAS surface observation reported 
thunderstorms, rainshowers, cumulonimbus, or lightning in either the significant weather 
section or the remarks section between 15 UTC and 22 UTC. These days were cross- 
referenced with the KSC/CCAS microburst-days identified by Sanger (1999), in his KSC 
microburst climatology, to determine which days had thunderstorms that produced 
microbursts and which days had thunderstorms that did not produce microbursts. 
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4. Methodology 


4.1 Introduction 

This thesis seeks to accomplish four objectives for the 45WS. 1) Fit a Weibull 
distribution to the 282 microburst gust speeds identified in Sangers’s (1999) KSC 
microburst climatology and integrate for the 45 WS convective wind warning criteria. 2) 
Evaluate MDPI using an expanded period of record. 3) Evaluate WINDEX as a predictor 
of microburst gust speed for the KSC/CCAS microburst forecast problem. 4) Develop a 
new microburst index using the MDPI and WINDEX parameters. The methodology 
presented in this chapter was designed to meet these objectives accurately and efficiently. 

4.2 Summary of Methodology 

1) Fit a Weibull Distribution to the KSC/CCAS observed microburst gusts 
identified by Sanger (1999) and integrate for convective wind warning criteria. 

2) Identify 1994 microburst cases and add to Sanger’s climatology. 

3) Distinguish microburst producing thunderstorm days from non-microburst 
producing thunderstorm days. 

4) Compute MDPI and WINDEX parameters. 

5) Evaluate MDPI using standard forecast verification techniques (i.e., 
probability of detection, false alarm rate, hit rate, threat score, Heidke skill score). 

6) Compare WINDEX versus observed microburst gust speed to determine 
predictive value. 
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7) Use WINDEX and MDPI parameters to develop a new microburst potential 


index by way of discriminant analysis. 

8) Compare new index to MDPI using standard forecast verification techniques. 


4.3 Fitting a Weibull Distribution to the KSC Microburst Climatology 

In 1939, the Swedish physicist Waloddi Weibull introduced a family of 
distributions now known as Weibull distributions. In some instances there are theoretical 
justifications for the use of a Weibull distribution, but in many applications they simply 
provide the best fit of observed data. The probability density function (pdf) for a Weibull 
distribution is 


f{x\y,a,/J) = \ 


- (x-r) a ~ ] e 

P“ 



0 


x>0 

x<0 


(?) 


where 

x is the range of the observed data 

/is the offset (equal to 0 if distribution begins at the origin) 
oris the shape parameter 
J3 is the scale parameter 

Weibull distributions are fit to a data set by adjusting the offset, shape and scale 
parameters to fit the distribution of the observed data (Devore 1995). Crystal Ball® was 
used to obtain the shape and scale parameters that fit the gust speed distribution of the 

(p\ 

282 microbursts identified by Sanger (1999). Mathcad was then used to integrate 


33 


equation (7) for the 45WS convective wind warning criteria yielding the probability of 
having microburst gusts greater than 35 knots, 50 knots, and 60 knots. Figure 7 is a 
sample Weibull distribution with y= 15, a= 2, and 30. 


Sample Weibull Distribution 



Figure 7. Sample Weibull Distribution with y = 15, a = 2, and p = 30. 


4.4 Identifying 1994 Microbursts Using WINDS Tower Data 

WINDS tower data was provided by the KSC for May - September 1994. A 
program written in IDL was used to identify microbursts based on Sanger’s (1999) six- 
condition algorithm since these microbursts would be added to his climatology. Three 
hundred seventy-two potential microbursts were identified using this approach. Further 
analysis of the data indicated that an overwhelming majority of these “microbursts” were 
a result of an apparent archival problem with the 1994 data. Thirty-four of the potential 
microbursts had peak gusts greater than 90 knots even though Sanger identified only 2 
microbursts with winds greater-than or equal-to 90 knots in his 4-year climatology. Most 
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of the extremely high peaks (70 knots or greater) had 5-minute mean winds, for the same 
period, of less than 10 knots. Additionally, many towers that have multiple sensors 
indicated a high peak wind at one level that was not reflected at levels above or below. 
Table 5 is a sample erroneous “microburst” from the 1994 WINDS data. Notice the peak 
wind of 117 knots at 54 feet with a mean during the same period of only 3 knots. Such a 
high peak wind should have a more significant influence on the mean for the same 
period. Additionally, the peak wind reading at 12 feet (just 42 feet below the 
“microburst”), for the same tower and period, is only 14 knots. 


Julian Date 

Time 

Tower 

Height 

Direction 

Mean 

Peak 

94142 

1215 

112 

12 

330 

3 

14 

94142 

1215 

112 

54 

333 

3 

117 


Table 5. Sample Erroneous Microburst. The 117-knot peak wind at 54 feet should have 
greater influence on the 5-minute mean for that height and the peak wind at 12 feet for 
the same period. 

“Microbursts” like the one described in Table 6 plagued the 1994 WINDS data, 
making it impossible to distinguish real microbursts from erroneous data. It is believed 
that these anomalies are an archival problem since the errors are not isolated to just a few 
towers and heights. As a result, the 1994 tower data was discarded and not added to 
Sanger’s (1999) climatology. 

4.5 Distinguishing Microburst/Non-Microburst Thunderstorm Days 

Surface observations from the KSC (KCOF) and the CCAS (KTTS), for the 
period May - September 1994-98, were used to identify thunderstorm days. A day was 
considered a “thunderstorm-day” if the KSC or CCAS surface observation reported 
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thunderstorms, rainshowers, cumulonimbus, or lightning in either the significant weather 
section or the remarks section between 15 UTC and 22 UTC. These days were cross- 
referenced with the microburst-days identified by Sanger (1999), in his KSC microburst 
climatology, to determine which days had thunderstorms that produced microbursts and 
which days had thunderstorms that did not produce microbursts. 

4.6 Computing MDPI and WINDEX Parameters 

The MDPI and WINDEX parameters were computed for microburst and non¬ 
microburst days that had 15 UTC soundings available. All parameters were computed 
using the equations and procedures described in Chapter 3. Of the 114 microburst-days 
identified by Sanger (1999), 92 had microbursts between 15 UTC and 22 UTC; of those, 
73 had 15 UTC soundings available. One hundred sixty-seven non-microburst producing 
thunderstorm days had thundertorms between 15 UTC and 22 UTC, and 15 UTC 
soundings available. Appendix E lists the MDPI and WINDEX parameters for the 240 
microburst and non-microburst producing thunderstorm days used in this study. 

4.7 Evaluating MDPI 

The MDPI was evaluated using the following forecast verification methods: 
probability of detection (POD), false alarm rate (FAR), hit rate (HR), threat score (TS), 
and Heidke skill score (HSS). These are the same methods used in previous MDPI 
studies conducted by the 45 WS. A chi-squared (x 2 ) test for goodness of fit was also 
computed for completeness. Additionally, a scatterplot of A0 e versus microburst 
occurrence was used to evaluate the MDPI critical threshold (CT). 
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Table 6 is a generic 2X2 contingency table used for forecast verification. Each 
position on the table displays the absolute frequencies, or counts, of each possible 
combination of forecast and observed pairs. Dividing these counts by the sample size 
(total number of forecast/observed pairs) transforms the counts into relative frequencies 
(Wilks 1995). The equations for computing the following accuracy and skill measures 
refer to Table 6. 

Observed 


Yes No 



Table 6. 2 X 2 Contingency Table for Forecast Verification (After Wilks 1995) 

4.7.1 Probability of Detection (POD) 

The probability of detection (POD), also known as prefigurance, is the fraction of 
occasions in which the forecast event occurred when it was forecast. For this study it is 
the likelihood that a microburst would be forecast, given that it occurred. The best 
possible POD is 1, indicating all microbursts that occurred were forecast. The worst 
POD is 0, indicating all microbursts that occurred were not forecast (Wilks 1995). In 
terms of Table 6 

POD = (8) 

a + c 
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4.7.2 False Alarm Rate (FAR) 


The false alarm rate (FAR) is the proportion of forecasts for which the forecasted 
phenomenon fails to materialize. For this thesis, it is the likelihood that a microburst is 
forecast but not observed. The best possible FAR is 0, indicating that all forecast 
microbursts were observed. The worst possible FAR is 1, indicating none of the forecast 
microbursts were observed (Wilks 1995). In terms of Table 6 


FAR = — 

a + b 


(9) 


4.7.3 Hit Rate (HR1 

The hit rate (HR), sometimes known as the proportion correct, is the fraction of 
the total forecasts that are correct. The worst hit rate is zero and the best hit rate is 1. 

The HR is sometimes multiplied by 100% and called the percent correct (Wilks 1995). In 
terms of Table 6 

hr = ^±± ( 10 ) 
n 

where n is the total number of forecast/observation pairs (a+b+c+d). 

4.7.4 Threat Score (TS! 

Sometimes it is not always desirable to use HR as a measure of forecast accuracy 
since it credits “y es ” and “no” forecasts equally. The threat score (TS), or critical success 
index (CSI), is often used when the event to be forecast, in this case microburst 
occurrence, occurs substantially less frequently than the non-occurrence. Essentially, the 
TS “is the number of correct “yes” forecasts divided by the total number of occasions on 
which that event was forecast and/or observed” (Wilks 1995). This is basically the same 
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as removing the correct “no” forecasts from the hit rate. The lowest attainable threat 
score is 0 and the best is 1 (Wilks 1995). In terms of Table 6 

TS = CSI = --—- (11) 

a + b + c 


4.7.5 Heidke Skill Score (HSS) 

Equations (8), (9), (10), and (11) measure the average correspondence between a 
forecast and the event the forecast is trying to predict. They are measurements of 
accuracy but not skill. Forecast skill is a measure of the relative accuracy of a set of 
forecasts, based on a set of control or reference forecasts. To have skill, the forecasts 
must perform better than the reference forecasts. The Heidke skill score (HSS) is a 
measure of forecast skill using the hit rate that would be achieved by random chance as 
the reference. It is the most commonly used skill score for summarizing 2X2 
contingency tables. Perfect forecasts receive a HSS of 1, forecasts that are as good as the 
reference receive a HSS of 0, and forecasts that perform worse than the reference receive 
a negative HSS (Wilks 1995). HSS is computed in terms of Table 6 by 


HSS = 


2 (ad - be) 

(a + c)(c + d) + {u + b)(b + d) 



To assess the validity of a forecast technique it is important to determine how the 


observed “hits” and “misses” in a 2 X 2 contingency table compare to the “hits” and 


“misses” that would be generated by a completely random selection process. The method 
for assessing the validity of the forecast technique is a chi-squared (x ) test for goodness 
of fit. The null hypothesis for the x 2 test is that the forecasts do not perform better than a 
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random selection process; the researcher’s hypothesis is that the forecasts do perform 

2 

better than a random selection process (Kachigan 1986). The equation for computing % 
is 


2 _ ^ (observed - expected) 
" aiiceiis expected 


where the expected value is computed by multiplying the marginal totals for the row and 
column, corresponding to the particular cell, and then dividing by the sum of all cells (n). 
In terms of Table 6 the expected value of cell a is 


expected value of cell a = 


(a + b)(a + c ) 
u + b + c + d 


If the x 2 value is greater than a critical value determined by an a priori significance level 
(a) then the null hypothesis is rejected and the forecast technique is said to be valid. If 
the x 2 value is less than the critical value the forecast technique is said to be without 
merit (Kachigan 1986). 

P-values were also computed based on the x 2 value. The /?-value is the area under 
the x 2 distribution to the right of the % 2 value. It is the smallest level of significance (a) 
at which the null hypothesis would be rejected. The p- value can also be thought of as the 
probability of getting a x 2 value greater than or equal to the one obtained by the sample 
data, given the null hypothesis is true. For example, if a x 2 value of 8.15 is computed 
based on a 2 X 2 contingency table, the associatedp-value of 0.0043 would indicate that 


there is a 0.43% probability of obtaining that value given the null is true (Devore 1995). 
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4.7.7 Evaluating MDPI Critical Threshold (CT) 


The MDPI critical threshold was evaluated using a scatterplot of A0 e versus 
microburst occurrence. The A0 e was plotted on the ordinate and microburst occurrence 
was plotted on the abscissa. Microburst occurrence was plotted as a 1 and non¬ 
occurrence was plotted as a 0. Figure 8 is an idealized scatterplot of A0 e versus 
microburst occurrence. Note that ideally there is very little overlap between the A0 e 
values corresponding to microburst occurrence and the A0 e values corresponding to non¬ 
microburst occurrence, with the overlap centered roughly on the MDPI critical threshold 
(3OK - horizontal line). 
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Figure 8. Idealized Plot of A0 e vs. Microburst Occurrence. Microburst occurrence is 
plotted as a 1 and non-microburst occurrence is plotted as a 0. Horizontal line marks the 
MDPI critical value (30K). 
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&v = £ (* -x)(yi-y) ( 16 ) 

i=i 

For the purpose of this thesis, 73 (n= 73) WINDEX/observed gust speed pairs 
{(xi,yd,~ • .,(x 73 ,y 7 j)} were evaluated. The WINDEX was computed using the 15 UTC 
soundings for microburst days identified by Sanger (1999). These WINDEX values were 
paired with the maximum microburst gust speed observed between 15 UTC and 22 UTC 
for the given microburst day. 

A scatterplot of observed maximum microburst gust speed versus WINDEX was 
used to visually determine the degree of linear correlation and to see if a non-linear 
relationship exists. Figure 9 shows an idealized scatterplot of observed maximum 
microburst gust speed versus WINDEX. The correlation coefficient for the idealized plot 
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is 1. Notice that the plot increases from left to right at a 45 degree angle showing a 
perfect linear relationship. 


Idealized Plot of 
Observed Microburst Gusts vs 
WINDEX 



WINDEX 


Figure 9. Idealized plot of observed maximum microburst gust speed versus WINDEX. 
The plot increases from left to right at a 45-degree angle indicating a perfect positive 
correlation (r = 1). 


4.9 New Index Development 

A statistical technique known as discriminant analysis was used to develop a new 


index for assessing wet-microburst potential for the Central Florida Atlantic Coast. 


Discriminant analysis is a method of classifying individuals or objects into mutually 
exclusive and exhaustive groups based on a set of independent predictor variables. It 


involves deriving linear combinations of the predictor variables that will discriminate 
between two or more previously defined groups such that misclassification error is 
minimized (Dillon and Goldstein 1984). This is done by assigning discriminant weights 
to the predictor variables. A discriminant score for each measured object is calculated by 
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multiplying the discriminant weight, associated with each predictor variable, by the 
object’s value on the predictor variable and then summing over the set of predictor 
variables. The equation for computing the discriminant scores is known as the 
discriminant function. For n number of objects to be classified and p predictor variables 
the discriminant function has the form 

Y=b'X (17) 

where 

Y is a 1 x n vector of discriminant scores 
b' is a 1 x p vector of discriminant weights, and 

X is ap x n matrix containing the values for each of the n objects on the p predictor 
variables. 

Objects are then grouped based on their discriminant scores and a discriminant 
cutoff score. Discriminant scores above the cutoff are placed into one group; 
discriminant scores below the cutoff are placed into the other. Figure 10 shows a 
graphical illustration of a two-group discriminant analysis with two predictor variables. 
The top part of the graph depicts a scatterplot of the two predictor variables (X! and X 2 ) 
for each object in the two groups. The lower part of the graph depicts the Y distributions 
(distributions of discriminant scores) for the two groups. The dotted line separating the 
two groups and two distributions is the cutoff score. Note that the greater the overlap of 
the two Y distributions (or scatterplots) the more errors in classification (Dillon and 
Goldstein 1984). 
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Figure 10. Graphical Illustration of Two-Group Discriminant Analysis. The upper part 
of the figure depicts a scatterplot (enclosed by ovals) of the two predictor variables (Xi 
and X 2 ). The lower portion of the figure depicts the distributions of discriminant scores 
(Y) for each group (After Dillon and Goldstein 1984). 


There were six steps to developing a new index using discriminant analysis. They 
were 1) choose predictor variables from the pool of MDPI and WINDEX parameters and 
other atmospheric measurements related to computing MDPI and WINDEX, 2) build a 
discriminant analysis input table (matrix) to train the discriminant model, 3) determine 
the discriminant weights, 4) compute discriminant scores for each group and determine 
the best cutoff score, 5) develop and evaluate the confusion matrix, and 6) compare the 
performance of the discriminant function (new index) to MDPI. 
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4.9.1 Choosing Predictor Variables 


For the purpose of this thesis, the two groups being discriminated were microburst 
producing thunderstorm days and non-microburst producing thunderstorm days. The 
predictor variables were drawn from the pool of MDPI and WINDEX parameters and 
other MDPI and WINDEX related variables. The nine MDPI and WINDEX parameters 
are described in Chapter 3. Other atmospheric parameters used as potential predictor 
variables were 

a. The pressure corresponding to the maximum 0 e in the lowest 150 mb of the 
atmosphere (mb) (P max fa) 

b. The pressure corresponding to the minimum 0 e between 650 mb and 500 mb 
(mb) ( P m in6k )• This is an indicator of the height of the cool, dry layer based on MDPI. 

c. The minimum 0 e regardless of the pressure level (K) (Min0 e . actu „i). This is the 
value used to compute A0 e in the MIST project (Atkins and Wakimoto 1991). 

d. The pressure corresponding to Min0 e . actua i (mb) (P min fa-actual)- This is an 
indicator of the height of the cool, dry layer based on the MIST project. 

e. The change in 0 e between the maximum 0 e in the lowest 150mb and Min0 e . ac tuai 
(K) (A0 e „ actua i)- This is the method used to compute A0 e in the MIST project (Atkins and 
Wakimoto 1991). 

Parameters for a stratified random sample of 50 microburst producing 
thunderstorm days and 50 non-microburst producing thunderstorm days were used to 
build (or train) the index. Equal samples of microburst and non-microburst producing 
thunderstorm days were used since the statistical software package used in this thesis (S- 
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Plus 2000) performs best using equal samples from each group. The remaining 23 
microburst producing thunderstorm days and a random sample of 39 non-microburst 
producing thunderstorm days were used to verify the new index. Since using 13 
parameters is cumbersome for an index, a method was developed to determine which 
parameters were most important for microburst development. 

A Pearson correlation coefficient (r) table (Appendix F), showing the r-values for 
each parameter as they relate to microburst occurrence and as they relate to each other, 
was developed using Statistix® for Windows®. The parameters that had the highest 
correlation to microburst occurrence and the least inter-correlation were used to build the 
index. They were A9 e . ac tmu Pminte-actmi, and 77 A discriminant analysis input table 
(matrix) was build using these parameters. 

4.9.2 Building the Discriminant Analysis Input Table 

A discriminant analysis input table is a tool for organizing the parameters that will 
be used to develop the discriminant function (Kachigan 1991). Table 7 is a generic 
discriminant analysis input table. Column 1 contains the objects that will be 
discriminated; for this thesis, it contains the days for which the predictor variables were 
measured. Column 2 indicates the group in which the object belongs. A 1 indicates a 
microburst producing thunderstorm day and a 0 indicates a non-microburst producing 
thunderstorm day. Columns 3, 4, and 5 are the predictor variables; they are equivalent to 
the transpose of matrix X in equation (17). 
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Object 

Group 

Predictor 1 

Predictor 2 

Predictor 3 


1 

Xu 

X12 

Xl 3 

Day 2 

1 

X21 

X22 

X23 

Day 3 

1 

X31 

X32 

X33 

Day 4 

1 

x 4 i 

X42 

X43 

Day 5 


X51 

X52 

X53 

Day 6 


*61 

X 62 

X 6 3 

Day 7 


x?i 

X72 

X73 

Day 8 


Xgl 

X 8 2 

X 83 

Day 9 

0 

x 9 , 

X92 

X 93 

Table 7. Generic Discriminant Analysis Input 1 

fable (After Kachigan 1991). 


Appendix G is the discriminant analysis input table for the 100 thunderstorm days 
used to train the discriminant function. Appendix H is the discriminant analysis input 
table for the 62 thunderstorm days used to verify the index. 


4.9.3 Determining the Discriminant Weights 

In 1936, Fisher developed a method for determining the discriminant weights 
based on the assumption that the variance-covariance matrices for each of the two groups 
have a common value (Dillon and Goldstein 1984). This common value is approximated 
by using the sample pooled variance-covariance matrix defined by the equation 

(».-l)[Sl] + (»2-l)[S»] 


Spooled — 


( Tl\ + Hi — 1 ) 


(18) 


where 

Spooled is the sample pooled variance-covariance matrix 

51 is the variance-covariance matrix for the predictor variables associated with group 1 

5 2 is the variance-covariance matrix for the predictor variables associated with group 2 
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nj is the number of objects measured in group 1, and 
n 2 is the number of objects measured in group 2 

The following equation is used to objectively determine the discriminant weights 
based on samples from each of the two groups. 

b = Spooled ' (Xi — X2) (19) 

where 

b is a vector of discriminant weights 

Spooled 1 is the inverse of the sample pooled variance-covariance matrix 
Xi is a vector of means of each predictor variable in group 1 
X 2 is a vector of means of each predictor variable in group 2 

S-Plus® 2000 was used to compute the sample pooled variance-covariance matrix 
and the group mean vectors for the stratified random sample of 100 thunderstorm days in 
Appendix G. A Mathcad template (Appendix I) was used to compute the discriminant 
weights based on the S-Plus® 2000 output. 

4.9.4 Computing Discriminant Scores and Determining the Cutoff 

A Mathcad template (Appendix I) was used to compute the discriminant scores 
for each group in the training sample in Appendix G using equation (17). The cutoff 
score was computed by the equation 

Ycutoff= b'^-^- (20) 

where b', Xi, and Xi are defined above. Equation (18) is valid only when the sample sizes 
of group 1 and group 2 are the same (Dillon and Goldstein 1984). Future observations of 
Y can be grouped based on Y cuio ff according to the rule 
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7 belongs to group 1 if it is greater than or equal to Y cuto ff 
Y belongs to group 2 if it is less than Y cuto ff 
4.9.5 Developing and Evaluating Confusion Matrix 

A confusion matrix is a type of contingency table used to evaluate a discriminant 
function’s ability to properly categorize objects or events. Table 8 is a generic confusion 
matrix for thunderstorm classification (microburst producing or non-microburst 
producing). Cells a and d indicate the number of correctly categorized thunderstorms. 
Cells b and c indicate the number of incorrectly categorized thunderstorms. Confusion 
matrices can be produced for the training data or for verification data (Dillon and 
Goldstein 1984). If the discriminant function is used to forecast events, the confusion 
matrix is the same as a contingency table for forecast verification. 


Actual Category 


o 

O) 

CD 
«#—* 
CO 

O 

c 

CO 

c 


o 

CO 



Table 8. Generic Confusion Matrix for Thunderstorm Category (Microburst Producing 
and Non-Microburst Producing). 


50 







4.9.6 Comparing Performance of the Discriminant Function (New Index) to MDPI 

Once the discriminant weights and the cutoff score have been computed equation 
(17) takes the form of a new index similar to MDPI. Future measures of the predictor 
variables can be used to compute a discriminant score (Y) that is categorized as either a 
microburst producing thunderstorm day or a non-microburst producing thunderstorm day 
based on the cutoff score (Y cu , 0 /f). Discriminant scores and MDPI values were computed 
using a stratified random sample of 62 thunderstorm days (23 that produced microbursts 
and 39 that did not) to compare the performance of the two indices. None of the 
validation thunderstorm days were used to train the discriminant function. Two 
contingency tables were produced: one containing the absolute frequencies, or counts, of 
forecast and observed microbursts based on the discriminant cutoff {Y cuto ff)\ the other 
containing the absolute frequencies based on MDPI. The indices were compared using 
POD, FAR, HR, TS, HSS and y\ Skill scores for POD, FAR, HR, and TS with MDPI as 
the reference forecast were also computed. Skill scores (SS) are interpreted as a 
percentage improvement over a reference forecast. They are computed with the 
following equation: 

SSref = A ' Aref xl00% (21) 

Aperf - Aref 

where 

SS re f is the skill score of a forecast technique with respect to some reference forecast 
technique 

A is the accuracy score for the set of forecasts associated with the forecast technique 
being evaluated 
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A re f is the accuracy score for a set of reference forecasts, and 

Aperf is the accuracy score that would be achieved by a set of perfect forecasts. 
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5. Results and Conclusions 


5.1 Introduction 

This chapter summarizes the results and conclusions of this thesis using the 
methods described in Chapter 4. The last section lists opportunities for future research. 

5.2 Fitting a Weibull Distribution to the KSC Microburst Climatology 

The program Crystal Ball was used to determine the Weibull shape and scale 
parameters that best fit the gust speed distribution of the 282 microbursts identified by 
Sanger (1999) in his KSC microburst climatology. The shape parameter {a) was 1.6263 
and the scale parameter {J3) was 16.88. Figure 11 shows a histogram of microburst gust 
speed with a rough distribution curve overlay. Figure 12 shows the fitted Weibull 
distribution based on the shape and scale parameters computed by Crystal Ball. Both 
figures were computed using Mathcad®. The curve in Figure 12 was integrated for each 
of the 45 WS convective wind warning critical thresholds yielding the probability of 
having winds greater than 35 knots, 50 knots, and 60 knots should a microburst occur. 

All computations were done using Mathcad®. The resulting probabilities are 

a. The probability of winds greater than 35 knots given microburst occurrence is .44 

b. The probability of winds greater than 50 knots given microburst occurrence is .08 

c. The probability of winds greater than 60 knots given microburst occurrence is .02 
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Figure 11. Histogram of Microburst Gust Speeds Identified by Sanger (1999) With 
Rough Distribution Curve Overlay. 



trace 1 


Figure 12. Weibull Distribution of 282 Microbursts Identified by Sanger. 
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5.3 Evaluation of MDPI 


MDPI was computed for 73 microburst producing thunderstorm days and 167 
non-microburst producing thunderstorm days using the procedure described in Chapter 3. 
Table 9 shows the 2 X 2 forecast verification contingency table for this data. The 
following accuracy and skill scores were computed using procedures described in 
Chapter 4. 

Probability of Detection (POD) = 65% Heidke Skill Score (HSS) = 8% 
False Alarm Rate (FAR) = 65% Chi-Square (x 2 ) = 2.12 

Hit Rate (HR) = 51% P-value = .145 

Threat Score (TS) = 29% 


MDPI Evaluation 


Yes 

c o 
CO 
O 

2 

o 

LL 

No 

Table 9. 2 X 2 Forecast Verification Contingency Table For MDPI 

Overall MDPI did not perform particularly well. A x 2 test was used to assess the 
validity of the index. The null hypothesis was that the “hits” and “misses” produced by 


Observed 

Yes No 


48 

93 

25 

74 
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MDPI are no better than those that would be assigned by a random selection process; the 
researcher’s hypothesis was that the “hits” and “misses” produced by MDPI are better 
than those that would be generated by a random selection process. A x 2 value of 2.706 or 
greater is required to reject the null for a significance level (a) of 0.1. The x 2 of 2.12 
indicates that the cell values produced by MDPI are not significantly better than those 
produced by random chance (the null is not rejected). Although the contingency table 
failed the % 2 test, the POD, FAR, HR, TS and HSS were computed since these are the 
measures used to evaluate MDPI in previous studies. The POD indicates that 65% of the 
microbursts would have been identified using MDPI, however, the high FAR indicates 
that 65% of the forecast “microbursts” would fail to materialize. This high FAR is 
unacceptable. The HR of 51% shows that using MDPI to determine if a microburst will 
occur is only slightly better than flipping a coin. The HSS, which compares MDPI to 
random guessing, indicates that MDPI only beats random guessing by 8%. 

The MDPI critical threshold (CT) was evaluated using a scatterplot of A0 e versus 
microburst occurrence (Figure 13). The large overlap of the two A0 e “populations” 
indicates that the current CT of 30K is as good as any other value depending on the level 
of acceptable error determined by the user. 
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A 8 e vs. Microburst Occurrence 



Figure 13. Scatterplot of A0 e vs Microburst Occurrence. A 0 indicates a non-microburst 
producing thunderstorm day and a 1 indicates a microburst producing thunderstorm day. 
The line marked CT is the MDPI critical threshold. 

5.4 Evaluation of WINDEX 

The Pearson correlation coefficient (r) was used to evaluate the degree of linear 
association between the maximum gust speeds predicted by WINDEX and the observed 
microburst gust speeds identified by Sanger (1999). Statistics® for Windows was used to 
compute r. The r -value was 0.24 (p-value 0.29) indicating very little correlation between 
WINDEX and the observed gust speeds. This result is consistent with r values computed 
in previous WINDEX studies (Roeder 1999). 

A scatterplot of observed maximum gust speed versus WINDEX (Figure 14) was 
used to visually determine the degree of linear correlation and to see if a non-linear 
relationship exists. It is clear from the plot that no discemable linear or non-linear 
relationship exists between WINDEX and the observed microburst gust speed. 
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Observed Gusts vs. WINDEX 



Figure 14. Scatterplot of Observed Microburst Gust Speed vs. WINDEX. It is clear 
from this plot that no discemable linear or non-linear relationship exists between 
WINDEX and the observed microburst gust speed. 

5.5 New Index for Assessing Wet-Microburst Potential 

A stratified sample of 50 microburst producing thunderstorm days and 50 non¬ 
microburst producing thunderstorm days was used to develop a new microburst potential 
index by way of discriminant analysis. The new index is 

„ T , b\(Maximum Oe- Minimum 0e) + bi(P min&) + 63 (E) ..... 

New Microburst Index =- (21) 

Ycutoff 

where 

b], b 2 , and hi are the discriminant weights 0.055 K' 1 , -0.007876 mb' 1 , and 0.966 km °C'’ 
respectively 

Ycutoff is the discriminant cutoff 3.231 

Maximum 6 e is the maximum value of 0 e in the lowest 150 mb of the sounding (K) 
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Minimum 9 e is the minimum value of 0 e in the entire sounding (K) 

PminSe is the pressure level corresponding to Minimum 6 e (mb), and 

Ais the environmental lapse rate from the surface to the height of the melting level (°C 

km 1 ). 

New index values of 1 or greater indicate a strong potential for microburst 
development should a thunderstorm occur. New index values of less than 1 indicate a 
low potential for microburst development should a thunderstorm occur. This follows the 
same convention as MDPI. 

Table 10 shows the confusion matrix for the 100 thunderstorm days used to train 
the index. Overall, the discriminant function classified 69% of the thunderstorm days 
correctly. The classification error rate was 31%. A x 2 test was performed on the 
confusion matrix to determine if the classification procedure, defined by the discriminant 
function and cutoff, is statistically significant based on the training data set. The % of 
14.49 and resulting /i-value of 0.0001 imply that the procedure is statistically significant. 


New Index Confusion Matrix 
Actual Category 



Microburst 

Non-Microburst 

it Category 

Microburst 

33 

14 

,i i 



i i 

O | 

CO 

3 1 

17 

36 


Table 10. New Index Confusion Matrix 
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A stratified random sample of 62 thunderstorm days (23 microburst producing 
and 39 non-microburst producing) was used to compare the performance of the new 
index to MDPI. None of the 62 verification thunderstorm days were used to build the 
new index. Table 11 shows the 2 X 2 forecast verification contingency table for MDPI 
based on the 62 thunderstorm day sample. 


^ Ves 

O') 

co 

o 

2 

o 

U_ 

No 


Table 11. 2X2 Forecast Verification Contingency Table for MDPI. This table is based 
on a stratified random sample of 62 thunderstorm days. 

The following accuracy and skill scores were computed from Table 11. 


POD = 69% 

HSS = 7% 

FAR = 60% 

tl 

H-k 

HR = 50% 

P-value = .522 

TS =34% 



Overall MDPI performed about the same as the evaluation presented in section 
5.3. Again a test was used to assess the validity of the index. The null hypothesis was 
that the “hits” and “misses” produced by MDPI are no better than those that would be 


MDPI Verification 
Observed 

Vies No 


16 

24 

7 

15 
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generated by a random selection process; the researcher’s hypothesis was that the “hits” 
and “misses” produced by MDPI are better than those that would be generated by a 
random selection process. Once again, the contingency table did not pass a x 2 test for a 
significance level of 0.1 suggesting that the MDPI shows little merit. The /7-value of .522 
suggests that the % 2 value of 0.14 would be attained 52.2% of the time as a result of a 
random selection process. Although the MDPI did not pass the x 2 test, POD, FAR, HR, 
TS, and HSS were computed for comparative purposes. The 2X2 forecast verification 
contingency table for the new index is presented in Table 12. 

New Index Verification 
Observed 



Yes 

No 

Vies 

-»—» 
c/) 

03 

O 

CD 

15 

16 

i_ 

O 



LL 



No 

8 

23 


Table 12. 2X2 Forecast Verification Contingency Table for the New Index. This table 
is based on a random sample of 62 thunderstorm days. 

The following accuracy and skill scores were computed from Table 12. 


POD = 65% 

HSS - 23% 

FAR = 51% 

X 2 = 3.39 

HR = 61% 

P-value = .066 

TS = 38% 
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2 

The new index outperformed MDPI in almost all accuracy and skill scores. A % 
test was used to assess the validity of the new index. The null hypothesis was that the 
“hits” and “misses” produced by new index are no better than those that would be 
generated by a random selection process; the researcher’s hypothesis was that the “hits” 
and “misses” produced by the new index are better than those that would be generated by 
a random selection process. The new index passed the % test for a significance level of 
0.1 (the null was rejected). The /(-value of 0.066 suggests that the x 2 value of 3.39 would 
be attained only 6.6% of the time as a result of a random selection process. This implies 
a significant improvement over the MDPI which reported a p-value of 0.522 for the same 
sample. 

As per equation 21, skill scores for the new index were computed using the POD, 
FAR, HR, and TS for MDPI as reference accuracy tests. The new index performed 13% 
worse than MDPI with respect to POD, 15% better than MDPI with respect to FAR, 22% 
better than MDPI with respect to HR, and 6% better than MDPI with respect to TS. The 
HSS shows that the new index beat random guessing by 23%. This is a marked 
improvement over MDPI which only beat random guessing by 7%. Table 13 compares 
the new index performance to that of MDPI. 
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MDPI 

New Index 

Skill Score 

POD 

69% 

65% 

-13% 

FAR 

60% 

51% 

15% 

HR 

50% 

61% 

22% 

TS 

34% 

38% 

6% 

HSS 

7% 

23% 

N/A 

X 2 

.41 

3.39 

N/A 

P-Value 

.522 

.066 

N/A 


Table 13. MDPI and New Index Performance Comparison. Skill scores were computed 
using MDPI as the reference forecast. 

5.6 Conclusion 

Overall neither the MDPI nor the WINDEX performed particularly well. Neither 
index should be used operationally at the KSC/CCAS. The MDPI showed very little 
improvement over random guessing (HSS = 7%) and did not pass a % 2 test for goodness 
of fit in either evaluation. The WINDEX showed very little correlation to observed 
maximum microburst gust speed (r = .24). Additionally, there was no discemable linear 
or non-linear relationship between WINDEX and observed maximum microburst gust 
speed. The Weibull distribution developed in this thesis is a good climatological starting 
point for assessing potential microburst gust speed. 

A new index for predicting microburst potential was developed by way of 
discriminant analysis. Although the new index outperformed MDPI, its performance 
was, by no means, stellar. The new index did pass a % 2 test for goodness of fit but 
showed only marginal improvement over random guessing (HSS = 23% and /7-value 
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.066). Further refinement of the new index is needed to make it a more useful operational 
forecasting tool. 

5.7 Opportunities for Future Research 

There are five recommendations for future reasearch. They are 

1) Add 1999 and 2000 microbursts to Sanger’s (1999) KSC microburst 
climatology and evaluate the new index using an expanded data set 

2) Redefine /’as the environmental lapse rate from the surface to the height of 
the minimum 0 e aloft and re-evaluate the new index 

3) Evaluate Tl, T2, and Snyder methods for forecasting microburst gust speed at 
KSC/CCAS (found in AFWA TN-98/002 Meteorological Techniques) 

4) Evaluate RAOB microburst prediction techniques used by the National 
Weather Service Office, Jackson, Mississippi for applicability at the KSC/CCAS. They 
include using Convective Available Potential Energy (CAPE), Vertical Totals (VT), and 
Wet Bulb Zero Height (WBZ) to assess microburst gust speed. 

5) Evaluate the possibility of downward advection of momentum as a possible 
cause for microbursts at the KSC. 
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Year 

Day 

Time 

Tower 

Height 

Direction 

Speed 

Gust 

1998 

3-Sep 

1815 

9001 

54 


14 

39 

1998 

3-Sep 

1825 

803 

12 

234 

9 

28 

1998 

3-Sep 

1830 

mm 

54 

244 

18 

40 

1998 


1835 

108 

MM 

245 

27 

43 

1998 


1900 

421 

54 

214 

27 

46 

1998 

7-Sep 

2145 

513 

30 

169 

19 

33 

1998 


1230 

311 

54 

67 

18 

29 

1998 


2055 

3 

54 

162 

25 

34 


71 
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Appendix B: Sample WINDS Data 


DAY TIME 


IDN 

Z 

DIR 

SPD 

GUS 

DDEV 

T 

TD 

RH 

94121 

0 

1 

6 





76 



94121 

0 

1 

12 

113 

6 

9 

16 




94121 

0 

1 

54 

123 

12 

13 

7 

76 



94121 

0 

3 

6 





78 



94121 

0 

3 

12 

99 

11 

13 

6 




94121 

0 

3 

54 

108 

13 

13 

7 

77 



94121 

0 

6 

6 








94121 

0 

6 

12 








94121 

0 

6 

54 








94121 

0 

6 

162 








94121 

0 

6 

204 








94121 

0 

108 

6 





75 



94121 

0 

108 

12 

100 

3 

7 

25 




94121 

0 

108 

54 

111 

8 

11 

11 

73 



94121 

0 

110 

6 





76 

70 

79 

94121 

0 

110 

12 

105 

5 

7 

11 




94121 

0 

110 

54 

111 

8 

9 

9 

76 

69 

79 

94121 

0 

110 

162 

101 

10 

11 

6 




94121 

0 

110 

204 

103 

11 

12 

5 

76 

69 

78 

94121 

0 

112 

6 





77 



94121 

0 

112 

12 

104 

4 

8 

15 




94121 

0 

112 

54 

102 

8 

10 

9 

77 



94121 

0 

303 

6 





76 



94121 

0 

303 

12 

108 

3 

4 

20 




94121 

0 

303 

54 

113 

6 

9 

13 

75 



94121 

0 

311 

6 





76 



94121 

0 

311 

12 

119 

5 

8 

16 




94121 

0 

311 

54 

105 

8 

10 

11 

74 



94121 

0 

313 

6 





76 

64 

68 

94121 

0 

313 

12 

127 

5 

8 

11 




94121 

0 

313 

54 

108 

7 

10 

10 

77 

66 

71 

94121 

0 

313 

162 

115 

11 

13 

6 




94121 

0 

313 

204 

115 

12 

13 

7 

76 

66 

73 

94121 

0 

313 

295 

120 

13 

14 

7 




94121 

0 

313 

394 

112 

13 

14 

7 




94121 

0 

313 

492 

109 

14 

14 

6 

75 

66 

73 

94121 

0 

403 

6 





76 



94121 

0 

403 

12 

134 

5 

7 

13 




94121 

0 

403 

54 

123 

6 

10 

11 

76 



94121 

0 

412 

6 





76 



94121 

0 

412 

12 

106 

5 

7 

14 




94121 

0 

412 

54 

98 

8 

11 

10 

75 



94121 

0 

415 

6 
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94121 

0 

415 

12 

124 

3 

5 

17 




94121 

0 

415 

54 

99 

6 

8 

10 

75 



94121 

0 

506 

6 





76 



94121 

0 

506 

12 

106 

4 

6 

14 




94121 

0 

506 

54 

107 

8 

10 

11 

76 



94121 

0 

509 

6 
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94121 

0 

509 

12 

103 

3 

4 

18 




94121 

0 

509 

54 

105 

5 

6 

9 

74 
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15 4 5 1998 1013 5 25.8 16.8 210 8 
15 4 5 1998 1011 22 24.6 17.6 999 999 
15 4 5 1998 1000 123 24.2 18.2 215 13 
15 4 5 1998 992 193 23.8 18.8 999 999 

15 4 5 1998 979.4 304 999 999 210 16 

15 4 5 1998 969 396 22 18.1 999 999 
15 4 5 1998 948 584 20.4 16.8 999 999 

15 4 5 1998 945.2 609 999 999 210 18 

15 4 5 1998 933 721 20.8 10.8 999 999 

15 4 5 1998 926 786 21.4 3.4 999 999 

15 4 5 1998 925 796 999 999 999 999 

15 4 5 1998 913 907 20.6 8.6 999 999 

15 4 5 1998 912.3 914 999 999 205 22 
15 4 5 1998 880.7 1219 999 999 205 22 

15 4 5 1998 850.2 1524 999 999 205 22 

15 4 5 1998 850 1526 14.8 4.8 210 22 

15 4 5 1998 837 1655 13.8 3.8 999 999 

15 4 5 1998 819.8 1828 999 999 210 20 

15 4 5 1998 790.3 2133 999 999 210 18 

15 4 5 1998 785 2189 10.6 -4.4 999 999 

15 4 5 1998 770 2349 9.8 -3.2 999 999 

15.4 5 1998 761.7 2438 999 999 215 18 

15 4 5 1998 733.8 2743 999 999 220 21 

15 4 5 1998 710 3013 4.6 -5.4 999 999 

15 4 5 1998 700 3138 3.8 -6.2 225 22 

15 4 5 1998 655.9 3657 999 999 240 22 
15 4 5 1998 632 3954 -3.5 -6.8 999 999 

15 4 5 1998 612 4207 -4.5 -6.7 999 999 

15 4 5 1998 607.3 4267 999 999 260 24 
15 4 5 1998 603 4323 -4.9 -8.7 999 999 

15 4 5 1998 590 4493 -6.3 -10.3 999 999 

15 4 5 1998 561.5 4876 999 999 260 26 
15 4 5 1998 550 5036 -10.7 -15.7 999 999 

15 4 5 1998 524 5405 -14.1 -16.5 999 999 

15 4 5 1998 518.6 5486 999 999 255 28 
15 4 5 1998 500 5770 -16.1 -17.5 250 28 

15 4 5 1998 490 5921 -17.5 -18.3 999 999 

15 4 5 1998 478.6 6096 999 999 245 29 

15 4 5 1998 463 6342 -19.9 -21.3 999 999 

15 4 5 1998 453 6503 -21.1 -23 999 999 

15 4 5 1998 442 6684 -22.1 -27.1 999 999 

15 4 5 1998 430 6885 -23.7 -31.7 999 999 

15 4 5 1998 419 7073 -25.3 -36.3 999 999 

15 4 5 1998 400 7410 -28.1 -38.1 265 34 
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Appendix D: IDL Program for Data Interpolation, and Computing MDPI and WIND EX 
pro parameters 10m, fii 

; This program will interpolate missing temp and dp data from upper-air soundings and 
; compute MDPI, WINDEX and associated parameters. Data is read-in in the format: 
;Time,Day,Month, Year,Press,Height,Temp,Dewpoint,Direction,Speed. All parameters 
;are saved to an output file. 

; Written By: Steven Dickerson 
; Last updated: 2 Dec 1999 

; Count the number of lines in the sounding 


n = 0 
s = " 
close, 5 
openr, 5, fn 

while not (eof(5)) do begin 
readf, 5, s 

if (strlen(s) GT 5) then n = n+1 

endwhile 

point lun, 5, 0 

; Now read in the sounding data 

data = fltarr(10,n) 
readf, 5, data 
close, 5 

;********This identifies the variable associated with each column of the array******* 

time = data[0,*] 
day = data[l,*] 
month = data[2,*] 
year = data[3,*] 
pres = data[4,*] 
hgt = data[5,*] 
temp = data[6,*] 
dp = data[7,*] 
dir = data[8,*] 
spd = data[9,*] 

;rh = data[10,*] 

; This identifies where the missing values of temp are 
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blanks = Where(strpos(temp, '999.0') GE 0, be) 
nonblank = Where(strpos(temp, '999.0') LT 0, nbc) 

for i = 0L, bc-1 do begin 

; find obs before 999.0 with number 

before = max(where(nonblank LT blanks(i))) 
after = min(where(nonblank GT blanks(i))) 

temp(blanks(i)) = temp(nonblank(before)) + ((temp(nonblank(after))- 
temp(nonblank(before))) * ((hgt[blanks(i)] - 

hgt[nonblank(before)])/(hgt[nonblank(after)]-hgt[nonblank(before)]))) 

;weight=( (hgt[blanks(i)] - hgt[nonblank(before)])/(hgt[nonblank(after)]- 
hgt[nonblank(before)])) 

;print,weight 

endfor 

; This identifies where the missing values of DP are 

blank = Where(strpos(dp, '999.0') GE 0, be) 
nonblanks = Where(strpos(dp, '999.0') LT 0, nbc) 

for i = 0L, bc-1 do begin 

; find obs before 999.0 with number 

before = max(where(nonblanks LT blank(i))) 
after = min(where(nonblanks GT blank(i))) 

dp(blank(i)) = dp(nonblanks(before)) + ((dp(nonblanks(after))-dp(nonblanks(before))) * 

((hgt[blank(i)] - hgt[nonblanks(before)])/(hgt[nonblanks(after)]-hgt[nonblanks(before)]) 

)) 

endfor 

. ************* This section builds a Theta-E profile of the Atms ***************** 
for i = 0, n-1 do begin 

; compute the temperature at the Lifted Condensation Level 
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TLCL = dp(0) - (0.212+.001571 *dp(0)-.000436*temp(0))*(temp(0)-dp(0))+273.16 

; compute the mixing ratio 

e = 6.112* exp(( 17.67 *dp)/(dp+243.5)) 
r = ,622*((e)/(pres-e)) 

; Now compute thetaE 

k=.2854 

thetaE = (temp+273.16)*((1000.0/pres) A (k*(l .0-(.28*r))))*exp(((3376/TLCL)- 
2.54)*r*(l+(.81*r))) 

endfor 

.************* -p^jg sec tion computes the delta-ThetaE needed for MDPI ********** 

position_lowest_150mb = where(pres GE (max(pres)-150)) 

lowPres = pres(position_lowest_150mb) 

maxThetaE = max(thetaE(position_lowest_150mb)) 

maxPres = lowPres(where(max(thetaE(position_lowest_150mb)))) 

;deltaThetaE=maxThetaE-min( thetaE) 

;position_minPres = where(thetaE EQ min(thetaE)) 

;minPres = pres(position_minPres) 

position minPres = where(pres GE 500 and pres LE 650) 
midlevelPres = pres(position_minPres) 
midlevel_thetaE = thetaE(positionminPres) 
minThetaE = min(midlevel_thetaE) 
minPres = pres(where(thetaE EQ min(midlevel_thetaE))) 

deltaThetaE = maxThetaE - minThetaE 

min_P_Total = pres(where(thetaE EQ min(thetaE))) 
min_TH_Total = min(thetaE) 
delta Total = maxThetaE - min TH Total 

. ***************** 'ptjjg section computes the Vertical Totals ****************** 
Pres850=where(pres EQ 850) 

Pres500=where(pres EQ 500) 


* 

;*Compute using min at mid* 
; * levels * 
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if ((min(Pres850)) AND (min(Pres500)) GT 0) then begin ; min statements needed to 
change from array to scalar 

VT = temp(where(pres EQ 850)) - temp(where(pres EQ 500)) 
endif else begin 

VT = 9999 ; if 850 or 500 mb level is missing 9999 is reported at VT 
endelse 

; ************** jhig section computes the WINDEX parameters ***************** 

.************** Find the Height of the Melting Level (Hm) in kilometers *********** 

position_below_zero = where(temp LE 0) 
hgt_below_zero = hgt(position_below_zero) 

Hm = hgt_below_zero(0)/1000 


.************* Find the Mix i n g Rat i 0 at t h e Melting Level (Qm) in g/kg *********** 

presbelowzero = pres(position_below_zero) 
pres_ML = pres_below_zero(0) 

dpbelowzero = dp(position_below_zero) 
dp_ML = dp_below_zero(0) 

vaporpress = 6.112*exp((17.67*dp_ML)/(dp_ML+243.5)) 

Qm = (,622*((vaporpress)/(pres_ML-vaporpress)))*1000 


;*** Find the Mean Mixing Ratio in the lowest 1km of the Atmosphere (QL) in g/kg *** 

position_lowest_lkm = where(hgt LE 1000) 
dplowestlkm = dp(position_lowest_lkm) 
preslowestlkm = pres(position_lowest_lkm) 


for i=0, n-1 do begin 

vaporpress = 6.112*exp((17.67*dp_lowest_lkm)/(dp_lowest_lkm+243.5)) 
mixingratio = (.622*((vaporpress)/(pres_lowest_lkm-vaporpress)))* 1000.0 

number = n_elements(mixingratio) 
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QL = (mixingratio(O) + mixingratio(number-l))/2 


endfor 




compute Rq 


jf* *JC jjC jjc j|* ^ 5f* #jc iji )j% )j{ ^ ?Jc rj{ ?|c 5|C ?JC 5j> 5|C 5j<S 5j> 5|C ^ ^ ^ ^ 


Rq = QL/12 


if (Rq GT 1) then Rq=1.0 


;***** Now find the Mean Environmental Lapse Rate (gamma) in degrees C per 1km 


temp_LT_zero=temp(position_below_zero) 
gamma = ((temp(0)-temp_LT_zero(0))/Hm) 

.********************* ]q ow Compute WINDEX (WI) ************************ 

WI = 5 * (Hm*Rq* (gamma A 2 - 30 +QL - 2*Qm)) A .5 

new_array = [time, day, month, year, pres, thetaE, hgt, temp, dp] 

output = [maxPres, maxThetaE, minPres, minThetaE, min P Total, min TH Total, 
deltaThetaE, delta_Total, VT, Hm, Qm, QL, Rq, gamma, WI] 
print, output 


openu, outfile, "parameters-micro-new.txt", /get_lun, /append 
printf, outfile, fn 
close, outfile 
freelun, outfile 

openu, outfile, "parameters-micro-new.txt", /get lun, /append, width = 300 
printf, outfile, output 
close, outfile 
free lun, outfile 

end 
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Appendix E: MDPI and WINDEX Parameters 
Legend 


YR -- Year 


MO — Month 
DY - Day 

OCR — Occurrence or non-occurrence of microburst. A 1 indicates that this day is a 
microburst producing thunderstorm day; a 0 indicates that this is a non-microburst 
producing thunderstorm day 

MaxPr — The pressure level (mb) corresponding to the maximum 0 e recorded in the 
lowest 150 mb of the sounding 

MaxThetE — The maximum 0 e (K) recorded in the lowest 150 mb of the sounding 

MinPr — The pressure level (mb) corresponding to the minimum 0 e found between 650 
and 500 mb 

MinThetaE — The minimum 0 e (K) found between 650 and 500 mb 

MinPActu — The pressure level (mb) corresponding to the minimum 0 e found at any level 
in the sounding 

MinTActu — The minimum 0 e (K) found at any level in the sounding 

MDPI — MaxThetE - MinThetaE 

DeltaActua -- MaxThetE - MinTActu 

Hm — The height of the melting level (km) 

Qm — The mixing ratio at the height of the melting level (g kg' 1 ) 

Q1 - The mixing ratio in the lowest 1 km of the atmosphere (g kg' 1 ) 

gamma — The environmental lapse (°C km' 1 ) rate from the surface to the height of the 
melting level 

WI — The computed WINDEX value 
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YR [ MO | DY | OCR | MaxPr l MaxThetE 1 MinPr 1 MinThetaE l MinPActu MinTActu MDPI DeltaActua Hm Qm Ql gamma Wl 
"95 5 11 1 ToTT 354.491 641 319.236 723 317.631 35.255 36.8599 5.364 1.0349 16.126 6.3386 57.006 

"95" ~5 12 0 1012 348.825 614 322.318 614 322.318 26.507 26.5071 4.875 1.7927 13.803 6.3348 49.799 

"95 5 i"9~ 0 1013 359.291 567 321,21 706 314.459 38.081 44.8315 4.875 0.6539 18.101 6.6765 61.831 

■95 5 23 1 Toi7 351.57 611 320.793 700 319.631 30.776 31.9383 4.267 2.2909 15.19 6.5756 50.436 



360.217 585 319.847 585 319.847 40.37 40.3701 4.673 0.5431 18.088 6.2059 54.596 

360.217 ~585 319.847 585 319.847 40.37 40.3701 4.673 0.5431 18,088 6.2059 54.596 

365.236 615 332.694 615 332.694 32.542 32.5415 5.178 5.2077 19.175 5.9869 43.477 


















































































































































YR | MO | DY | OCR MaxPrI MaxThetE I MinPr I MinThetaE MinPActu MinTActu MDPI DeltaActua Hm Qm Ql gamma Wl 
"95 8 23~ 1 1011 362.096 545 337.683 545 337.683 24.413 24.4132 5.193 5.6465 17.653 5.5844 31.3 

"95 8 26 1 1013 361.419 543 335.866 463 332.586 25.553 28.8326 4.87 5.8155 17.231 5.9548 36.69; 

"95 8 27 1 ~1014 353.264 603 ~ 325.745 603 325.745 27.52 27.5195 5.18 1.9308 14.828 5.7874 43.272 

”95 8 30 0 Tori 367.383 500 327.931 500 327.931 39.452 39.4521 5.18 3.599 20.104 5.6854 44.411 
















































































































































































YR | MO DY OCR MaxPr MaxThetE MinPr MinThetaE MinPActu MinTActu MDPI DeltaActua Hm Qm _ Ql gamma Wl 

"96 6 23 0 1016 361.159 598 323.041 598 323.041 38.1 IB 38.1181 5.42 0.4481 17.37 6.08B6 56.4B4 

"96 6 25 0 1016 369.48 589 320.546 589 320.546 48.9 34 48.9344 4.877 0.674 18.7 6.4678 59.652 

"gg 6 26 i 1017 358.638 621 323.643 621 323.643 34.995 34.9947 4.414 3.2057 16.378 7.0229 56.851 

"96 6 28 0 1015 358.811 572 333.673 661 333.348 25.138 25.4626 4.806 4.6229 17.799 5.826 38.746 














































































































































































































































































































YR | MO DY OCR MaxPr MaxThetE MinPr MinThetaE MinPActu MinTActu MDPI DeltaActua Hm Qm Ql gamma 
"97 6 28 0 1018 356.278 647 325.338 714 313.413 30.939 42.8651 4.543 3.9899 15.66 6.4518 

~97~ 6 ~29 1 TtmT 350.773 500 313.487 500 313.487 37.286 37.2857 4.442 3.7204 14.901 6.5961 

97 6 4 0 1007 352.403 643 319.711 715 316.165 32.691 36.2374 4.703 1.4529 14.25 5.8898 

"97 6 5 0~ 1007 351.35 643 317.413 785 312.292 33.936 39.0573 4.535 0.5422 13.246 6.0858 






















































































































































YR | MO | DY | OCR 1 MaxPrl MaxThetE I MinPr I MinThetaE MinPActu MinTActu MDPI DeltaActua Hm Qm Ql gamma Wl 

"97 7 30 0 1019 363.898 542 329.290 689 328.662 34.599 35.2354 4.713 4.8842 16.797 6.5774 48.896 

97 7 ~31~ 0 1019 356.154~ 579 325.051 725 323.502 31.104 32.6526 5.231 1.5411 16.678 6.4319 57.138 

"97 8 IfT 0 1022 356.967 618 324.611 618 324.611 32.355 32.3553 5.041 2.9512 16.036 6.1495 47.562 

"97 8 n 0 1022 358.513 580 324.525 488 322.35 33.988 36.1621 4.78 1.8974 17.136 5.8579 45.933 









































































































































































YR MO DY OCR MaxPr MaxThetE MinPr MinThetaE MinPActu MinTActu MDPI DeltaActua Hm Qm _ Ql gamma Wl 

~98~ 5 1~9~ 0 ~ 017 346.26 643 315.082 865 ~311.251 31.177 35.0082 5.142 0.3371 14.42 6.1439 52.563 

98 5 ~28~ 0 1017 358.478 603 335.031 850 324.551 23.447 33.9265 5.152 5.5947 13.948 5.8229 29.299 

98 6 18 0 1018 348.704 643 322.921 700 321.317 25.783 27.3869 4.951 1.6463 13.936 6.2615 49.567 

"98 6 19 0 - 1017 355.876 572 322.318 572 322.318 33.559 33.5588 4.571 2.3196 13.151 7.4997 63.024 
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Appendix F: Pearson Correlation Coefficient (r) Table 


STUDENT EDITION OF STATISTIX 


CORRELATIONS (PEARSON) 


HM 

OCCUR 

-0.0583 

HM 

MDPI 

MDPI 

0.2385 

0.0415 


DELTAACTU 

0.2636 

0.1336 

0.8505 

GAMMA 

0.2010 

-0.5052 

0.3540 

MINPACTUA 

-0.3196 

-0.0833 

-0.0701 

MINPRES 

-0.3418 

-0.0942 

-0.2237 

QL 

0.0548 

0.4661 

0.4072 

QM 

0.1321 

-0.0511 

-0.2159 

WI 

0.0899 

-0.0698 

0.5863 

MAXTHETAE 

0.1846 

0.5080 

0.6224 

MINTHACTU 

-0.1263 

0.3781 

-0.3786 

MINTHETAE 

-0.0944 

0.4888 

-0.5563 

QM 

QL 

0.2745 

QM 

WI 

WI 

0.0330 

-0.5209 


MAXTHETAE 

0.8017 

0.2058 

0.2488 

MINTHACTU 

0.3889 

0.4122 

-0.3799 

MINTHETAE 

0.3555 

0.4813 

-0.4495 


CASES INCLUDED 100 MISSING CASES 0 



12/21/99, 2:04:22 PM 

DELTAACTU GAMMA 

MINPACTUA MINPRES 


0.3001 




-0.0996 

0.0824 



-0.2134 

-0.0234 

0.6399 


0.3803 

-0.2189 

-0.1603 

-0.2120 

-0.1572 

-0.0662 

-0.3535 

-0.3106 

0.5264 

0.7753 

0.2067 

0.0418 

0.6272 

0.0075 

-0.1971 

-0.2233 

-0.5568 

-0.3598 

-0.0882 

0.0234 

-0.3693 

-0.4229 

-0.1239 

0.0352 


MAXTHETAE MINTHACTU 


0.2977 

0.3041 0.7770 
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Appendix G: Discriminant Analysis Input Table for Training Data 


Number 

Occur 

DeltaActua 

MinPActua 

Gamma 

1 

1 

30.5078 

535 

6.13497 

2 

1 

30.6852 

614 

6.79737 

4 

1 

35.6141 

627 

6.74633 

6 

1 

25.5138 

477 

5.44959 

8 

1 


545 

5.58444 

10 

1 

34.9978 

628 

6.84346 

11 

1 

42.0995 

561 

6.54879 

12 

1 

28.8326 

463 

5.95483 

13 

1 

32.9438 

705 

6.56406 

14 

1 

27.5195 

603 

5.78737 | 

15 

1 

34.4821 

700 

6.26058 

18 

1 




19 

1 


678 

6.27436 

20 

1 

31.7443 

541 


22 

1 

29.1636 

666 

5.99342 

25 

1 

27.7669 

700 

6.63658 

27 

1 

33.1595 

602 

6.05701 

31 

1 

41.9129 

688 

6.56737 

32 

1 

44.5963 

705 

6.91173 

33 

1 

36.4901 

521 

6.24724 

34 

1 


555 

5.99642 

35 

1 


593 

6.55659 

36 

1 

42.9333 

400 

6.15171 

37 

1 

37.7194 


6.56514 

38 

1 

12.8243 

631.9 

5.17449 

39 

1 

45.0815 

400 

6.23343 

40 

1 

36.6702 

622.9 

6.25134 

41 

1 

32.7751 

629.8 

7.01451 

42 

1 


500 

6.35015 

43 

1 


500 

6.94317 

44 

1 

28.1281 

579.8 

6.72157 

48 

1 

—klllifeM 

400 

5.71873 

49 

1 



6.94082 

50 

1 



6.40889 

51 

1 

37.2857 

500 

6.59613 

52 

1 

47.4426 

500 

6.57978 

53 

1 

43.9083 

400 

6.07468 

54 

1 



6.49458 

56 

1 



6.76277 

57 

1 


587.8 

6.05782 

58 

1 

58.4116 

500 

6.10432 

59 

1 

50.0744 

500 

6.20872 

60 

1 

35.8315 

621 

6.32357 
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121 

o 

1 29.2216 

705 


90 


Gamma 



BaEtoai:« 

6.68118 

6.91585 

7.666 

6.15315 

6.89302 

6.09223 

6.55219 

6.30631 

5.98687 

6.1672 

6.43821 

5.66572 

5.93245 

5.68543 

6.85841 

6.3762 

6.21611 

5.89278 

7.0373 

6.01418 

6.50185 

6.90954 

5.96013 

6.12171 

6.60046 

6.40298 | 

5.87224 

6.0095 

5.54828 

6.06353 

6.46421 

6.62836 

6.82864 

6.88285 

6.45184 

6.14228 



5.91333 

6.39273 

6.82418 

6.14947 

5.93137 


6.29795 

























































































































Number 

Occur 

DeltaActua 

MinPActua 

Gamma 

124 

0 


451 


125 

0 

32.3209 

650 


131 

0 

32.6526 

725 

6.43193 

138 

0 

37.9366 

669 

6.14901 

142 

0 

31.1764 

649 


144 

0 


654 


145 

0 


567 


152 

0 

20.0883 

423 

5.52596 

159 

0 

31.231 

642 

5.96217 

161 

0 


637 

6.69228 

163 

0 


700 

6.0552 
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Appendix H: Discriminant Analysis Input Table for Verification 


DeltaActu 


36.8599 

34.2638 


29.8398 


31.9383 


MinPActu 


723 6.33855 

571 6.38257 


697 7.10294 

700 6.57558 


41.0621 


23.2509 


36.7213 


48.1011 


41.2154 


29.3167 


32.8301 

38.4923 


25.2333 

28.3034 


37.6819 


33.1985 


38.1181 


30.7538 


42.43 


18.0734 


32.6949 


37.7988 


34.0359 


37.571 


36.1621 





6.77486 


5.58692 




619.5 

6.32353 

500 

5.93617 



2 

6.56379 

0 

6.42412 



6.97294 


6.75672 



6.46495 




5.91607 


5.90926 


5.81623 


6.29317 

6.98701 




6.26208 


6.08861 


6.68811 




6.57792 


6.39768 


5.90724 


6.50294 


5.8579 
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143 

0 


716 

5.78753 

146 

0 


639 

5.87884 

154 

0 

41.5662 

591 





614 


162 

0 


571 


164 

0 

30.4287 

656 

5.86334 

165 

0 

27.4828 



166 

0 

24.8502 

439 

5.63372 
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Appendix I: Discriminant Analysis Mathcad Template 


ORIGIN e1 

This is a template that uses output from S-Plus 2000 to determine the cut-off between 
microburst occurrence and non-microburst occurrence. 


36.29817 


32.17186 


xbar j := 371.37 xba^ : = 632.90 


These are the group means. 


6.330973 


6.2083707 


38.16636 -11.764 .847691 

Spooled := - n -764 8488.621 6.187427 Pooled Sample Variance Covariance Matrix 

[.847691 6.187427 .180611 

b hat :=s pooled' 1 ( xbaj: 1 - xbar 2) This C0IT1 P utes discriminant vector (weights). 


b hat = -7.876 10’ 
0.966 


s the discriminant vector (sample approximations of discriminant weights) 


Cutoff :=b i 


( xbar i + xbar 2 ) This computes the cuttoff 


Cutoff = (3.231) 


A:\verify62.txt 

DeltaActual MinPActual V :=x 

Y : = b h a t -DeltaActual + b ^ at -MinPActual + b 


r This is the discriminant 
function 



1 

1 

2.466 

2 

3.563 

3 

3.022 

4 

2.604 

|5 : i 

3.827 

6 

3.619 

ff- 

4.166 

8 

3.828 

fs. 

3.333 

io 

1.642 

11 

4.045 

12 

3.743 

13 

4.621 

14 

2.658 

15 

3.503 

16 

2.56 


Cutoff = (3.231) 

The first 23 are microburst days the second 39 are non-microburst days 
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